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Abstract 


With mass uptake of mobile broadband such as HSPA, LTE and advanced LTE there 
will be a rapid increase of demand for capacity in the mobile backhaul. The traditional 
frequency bands used for mobile backhaul microwave links are narrow and may become a 
future capacity bottleneck. Line-of-sight MIMO is one possible technology for increasing 
the spectral efficiency of these links. By using multiple antennas at both the transmitter 
and receiver the capacity of such systems could be multiplied using the same spectrum. 


Line-of-sight MIMO works by making the signal streams orthogonal to each other 
either by spatially separating the antennas at the transmitter and receiver or by using 
orthogonally polarized antennas. 


In this thesis both 2x2 spatially-separated and 4x4 spatially-separated and dual-polarized 
line-of-sight MIMO are numerically investigated. Adaptive algorithms for a space-time- 
equalizer are proposed and evaluated against detrimental affects such as frequency selective 
channels, mast swing (antenna movement), polarization leakage and phase noise (corre- 
lated and differential). The algorithms are either decision-directed or blind, those are 
LMS and WRLS and CMA, SCMA and MCMA respectively. 


The simulations show that line-of-sight MIMO meeds higher system margins and a 
more complex receiver structure. Convergence of blind algorithms can be difficult to 
achieve for some channel conditions as well as for differential phase noise. Mast swing 
will further more put higher requirements on spatially-separated MIMO systems than on 
its SISO counterpart. 
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Chapter 1 


Introduction 


1.1 Background 


Todays mobile traffic is dominated by voice from GSM and WCDMA. In the not-so distant 
future the traffic will be dominated by broadband data. The evolution of HSPA and the 
arrival of LTE will increase data rates of future wireless networks manyfold creating tougher 
requirements on the capacity of the backhaul networks. A backhaul link is a link that connects 
base stations to network controllers, these base stations are located on the field often making 
them inaccessible to fiber- and DSL solutions. A microwave radio link is thus a commonly 
used solution since it does away with a direct cable connection. It is very important that the 
throughput of these radio links keep up with demand, preferably without requiring a purchase 
of more spectrum which can be scarce and expensive. 

Wireless MIMO technology allows for higher throughput for the same spectrum because of it's 
ability to transmit and receive multiple streams simultaneously. For most applications MIMO 
requires a non-line-of-sight- and a rich scattering environment to work properly but microwave 
links usually have non of this. These links use highly directional antennas that have a direct 
line-of-sight (LoS) component dominating the transmission. However, it may be shown that 
at microwave frequencies it is possible to transmit multiple streams over a single frequency 
band if the antennas at the transmitter and receiver are separated sufficiently. Combining this 
with dual-polarized antennas it is possible to quadruple the throughput without making the 
required antenna separation infeasible. Unlike the mobile MIMO channel which is inherently 
random the LoS MIMO channel is determined by geometry which essentially remains static 
over the operational lifetime. A properly installed link is thus more predictable than its 
mobile counterpart. 


1.2 Problem Description 


This thesis was started to be able to better understand what requirements the LoS MIMO 
channel would put on the system, both spatially separated and dual polarized systems. Adap- 
tive algorithms for these systems were investigated, e.g. equalizers able to cancel inter-stream 
interference and cross polarization distortion as well as conventional intersymbol interference. 
The performance of these algorithms were investigated for different channel conditions having 
time-variant frequency selectivity due to mast swing. Phase noise tolerances were also inves- 
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tigated as well as the effects of suboptimal antenna installation and high cross polarization 
couplings. The algorithms considered were training-based, semi-blind and blind. The thesis 
outcome should give a guideline of tolerances of such systems to channel effects as well as on 
installation requirements. 


1.3 Thesis Outline 


The report consists of seven chapters and two appendices. The first chapter along with the 
abstract give an overview of the work carried out. The second chapter is a high level overview 
of wireless communication using microwaves. Chapter three discusses MIMO technology with 
an emphases on LoS MIMO both spatially-separated and dual-polarized. Chapter four is a 
theoretical overview of equalizers and adaptive algorithms for MIMO. Chapter five describes 
how the simulations were carried out. Chapter six lists all the results obtained by simulations 
and finally chapter nine summarizes the results and gives some conclusions. 


1.4 Notation 


1.4.1 Abbreviations 


ADC Analog to Digital Converter 
ADSL Asymmetric Digital Subscriber Line 
AWGN Additive White Gaussian Noise 


BER Bit Error Rate 

CIR Channel Impulse Response 

CMA Constant Modulus Algorithm 

CSIR Channel State Information at the Receiver 
CSIT Channel State Information at the Transmitter 
DAC Digital to Analog Converter 

DD Decision Directed 

DFE Decision Feedback Equalizer 

DP Dual Polarized 

DTFT Discrete-time Fourier transform 

FIR Finite Impulse Response 

GUI Graphical User Interface 

IEEE The Institute of Electrical and Electronics Engineers 
ISI Intersymbol Interference 

LMS Least Mean Square 

LoS Line of Sight 

LS Least Square 

LTE Long Term Evolution 


MCMA Modified Constant Modulus Algorithm 
MIMO Multiple Input Multiple Output 

MISO Multiple Input Single Output 

MLSE Maximum Likelihood Sequence Estimation 
MMA Modified Modulus Algorithm 

MRLA Modified Recursive Least Square 
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MSE 
MTI 
NLoS 
PRBS 
PSK 
QAM 
RLS 
Rx 
SCMA 
SIMO 
SIP 
SISO 
SNR 
SP 
STE 
SVD 
Tx 
WiMAX 
WRLS 
XPD 
XPIC 


Mean Square Error 

Multiple- Transmitting-Interference 
None-Line of Sight 

Pseudorandom Binary Sequence 

Phase Shift Keying 

Quadrature Amplitude Modulation 
Recursive Least Square 

Receiver 

Simplified Constant Modulus Algorithm 
Single Input Multiple Output 
Superimposed Pilots 

Single Input Single Output 

Signal To Noise Ration 

Single Polarized 

Space Time Equalizer 

Singular Value Decomposition 
Transmitter 

Worldwide Interoperability for Microwave Access 
Weighted Recursive Least Squares 
Cross-Polarization Discrimination 
Cross-Polarization Interference Canceler 


1.4.2 Symbols 


d SV ooo 
y 
* 


The mixing parameter for SCMA 
Roll-off-factor for the root-raised Cosine filters 
Maximum tilt of antenna masts 
System delay in symbol periods 
Estimated Signal (Hard Decision) 
Condition Number 

Singular Value 

The ground reflection factor 

The stepsize of the CMA. 

The stepsize of the LMS in the Decision-Directed Mode 
The step-size of the MCMA 
Angular frequency 

Carrier Angular Frequency 
MIMO Channel Matrix 
Coherence bandwidth 

Signal bandwidth 

Capacity 

Eigenvalues 

Rx antenna separation 

Tx antenna separation 

Frequency of mast swing 
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Fe Carrier frequency 

F; The sampling frequency 

k Wavenumber 

À Weighting Factor for WRLS 

L Number of filter taps for each filter in the STE 
Li 


r Number of training symbols 
Ly Tx and Rx antenna mast horizontal distance 
Lr, Rx landscape height 
Ln Rx antenna mast height 
Li The placement of the central tap in each filter in the STE 
Lr Tx antenna mast height 
M Number of Rx antennas 
N Number of Tx antennas 


p(t) Matched Filter 

P, Bit Error Rate 

Q Oversampling rate at the Tx filter 

R The average constant magnitude of Blind algorithm 
Ry Rank of matrix H 


Tn Received Symbol n 

Sn Transmitted Symbol n 

Tos Symbol period of the oversampled signal 
T Coherence time 

T; Symbol period 

VMs Velocity of the antennas 

Yn Equalized Signal (Soft Decision) 

M Modulation order for M-QAM 


1.4.3 Operators 


eu Transpose 
(9H Complex conjugate transpose (Hermitian transpose) 
R{-} Real-part operator 
S{-} Imaginary-part operator 
l 


E[-] Expectation value operator 
trace{:} The trace of a matrix 
«5:2» Inner product 

EC) Laplace transform operator 
ó(-) Dirac's delta function 
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Chapter 2 


Wireless Microwave Systems 


This chapter will give a high level overview of wireless digital communication. Basic system 
designs are discussed as well as challenges that need to be overcome by these systems. 


2.1 Overview 


Communication devices have become so customary and their use so habitual that we hardly 
ever take time to wonder how all of this can even be possible. The complexity of todays sys- 
tems are immense but their basic design and operational principles have remained the same 
for quite some time. This design and these principles can easily be understood at a high level. 


Be it wireless cell phones or wired ADSL computer modems the purpose is the same, to com- 
municate through an exacerbating channel. In order to do that the transmitter sends special 
signals that are tailored in a way so that they are retrievable at the receiver. The receiver 
may improve the quality of the received signal by applying signal processing in order to filter 
out noise, remove interfering signals and combat other depriving effects of the channel to get 
an estimate of the symbols being sent. 


Analog and digital communication differ in the fact that digital signals have a finite signal 
set while analog transmission has an infinite amount (is continuous). Wherever one desires 
digital information at the endpoints of the system it is in fact necessary to use digital com- 
munication. Seeing that the vast majority of information is stored and processed digitally 
and the fact that digital communication can be designed to be functionally superior in almost 
all cases compared to analog communication, it is no wonder that the vast majority of to- 
day’s communication systems are digital. From now on only digital systems will be considered. 


For all communication systems there are two universal design parameters that need to be 
optimized. These are spectral efficiency and power efficiency, both have a profound impact 
on wireless links. Spectral efficiency quantifies the information content transmitted over a 
given bandwidth and has the unit bits/Hz/s. Power efficiency quantifies the signal-to-noise 
(SNR) requirement at the receiver to obtain a specific bit error rate (BER). Design and op- 
timization needs then be done without adding excessive complexity to the system. One can 
state that the brunt of the research effort done in recent times on digital wireless transmission 
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is concerned with optimizing for these three requirements. 


A high level view of a communication system can be stated as follows. The transmitter takes 
data and maps it into a sequence of suitable modulated and filtered symbols that are sent 
over a channel. The channel, which is comprised of a certain environment, then distorts the 
signal by adding noise and interfering signals, spreading out the signal in time causing inter 
symbol interference (ISI) and attenuating the signal power. The receiver then needs to undo 
as much of the effects of the channel as possible before interpreting the transmitted data and 
outputting (hopefully) the original data. 


2.2 Basic System Design 


A schematic representation of a communication system can be seen in Figure 2.1. Each block 
represents a conceptual or logical part of the system and is useful in obtaining a high level 
understanding. 


Source 


Source/ 
Channel Modulator 
Encoder 


Channel 


Source/ 
Channel Demodulator 
Decoder 


Figure 2.1: A schematic view of a basic communication system. 


The transmitter is fed by some signal source (e.g. a voice signal) and the source encoder 
converts that signal into a suitable digital format. The channel encoder then adds redun- 
dancy (i.e. extra bits) to the bit sequence to enable error detection and/or error correction at 
the receiver. The modulator then maps the bit sequence into discrete symbols. There are a 
number of possible symbol sets but two of the more popular can be seen in Figure 2.2(a) and 
Figure 2.2(b). The symbols (S;) are mapped onto a 2-D symbol space with an orthogonal 
basis expressed using real- and imaginary numbers. Having the symbols in complex format 
is a convenient expression, in a real transmitter the imaginary- and real parts are separated 
before being signal processed and finally mixed up to carrier frequency using the same carrier 
but with a 90? offset (in phase and quadrature phase). This phase offset in the carrier realizes 
the orthogonal basis expressed by the complex numbers. Transmission scheme like this are 
called quadrature multiplexed and part of such transmitter can be seen in Figure 2.2(c). The 
symbol S; is split into its real and imaginary parts before being processed and mixed up in 
frequency using carriers with a 90? offset. 


The pulse shaping filter is denoted as p(t), it maps the symbols onto a sequence of Nyquist 
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pulses. If the time between symbols is denoted as 7; (symbol time) then a Nyquist pulse 
must satisfy for some integer N the relation 


1 ifN=0 
pT.) = | 0 if N40 


This means that the pulse only has a value at the current sampling time but is zero for all 
other sampling times. This can be seen in Figure 2.2(d) where three symbols have been sent 
through a Nyquist pulse shaping filter. At exactly the symbol time (multiples of the symbol 
period) the surrounding symbols are zero so no interference between symbols is observed. 
After the signal has gone through the pulse shaping filter at the transmitter, and turned 
analog using a DAC (Digital to Analog Converter), it is up converted using a mixer. 


A mixer is an essential part for most wireless systems. A typical digital radio performs 
its signal processing at the baseband, that is between 0 and B [rad/s]. In order to utilize 
the available spectra in an efficient way, the signal is upconverted to a frequency band pre- 
allocated for the given application. It is also impractical to transmit a baseband signal so it 
is necessary to upconverted to the carrier frequency we. Since all signals can be expressed 
via Fourier analysis by a sum of sinusoids it is enough to consider only a signal z(t) = 
Ay sin(w»t) where 0 € w < B and Aj, is the amplitude. Mathematically a mixer performs 
the upconversion in the following way 


Ze (t) = z(t) sin (wet) = A [cos((we — w)t) — cos( (we + wp)t)] (2.1) 


where sin(w-t) is generated by an oscillator. From the expression in equation (2.1) it can be 
seen that the signal frequency is now spaced between we — B and we + B. 


Imt) Im{} 


Ref} 


b) d) 


Figure 2.2: a) 16-PSK modulation format (signal constellation) where S; represents one 
symbol b) 16-QAM modulation format c) A quadrature multiplexed transmitter. d) 3 symbols 
that have gone through a pulse shaping filter. 


At the receiver the signal is down converted with a mixer, filtered with a matched filter (same 
kind of filter as at the transmitter) to minimize the effect of noise and then sampled with 
an analog to digital converter (ADC). Then it is sent through the demodulator followed by 
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source/channel decoders which undoes the encoding operations at the transmitter. Finally 
the data sink operates on the data in some way, e.g. reproduces a voice signal. A much more 
detailed discussion can be found in many textbooks such as [1, 53]. 


2.3 The Wireless Channel 


A channel is here defined as the paths over which an electrical signal traverses between a 
transmitter (Tx) and receiver (Rx). That can include parts of the Tx and Rx, such as am- 
plifiers and oscillators, depending on preference. Wireless channels are specially challenging 
compared to wired channels. The ether (or transmission media) is firstly shared by all users 
making it very prone to interference, secondly it is composed of the ever changing environ- 
ment we all live in making it highly time varying. From a signal processing stand point the 
channel can introduce noise and other interfering signals as well as inter symbol interference 
(ISI) and fading. These effects will be discussed in the following sections. 


2.3.1 Noise 


Additive White Gaussian Noise (AWGN) is a common noise model for numerical and ana- 
lytical analysis since it approximates many noise sources (e.g. thermal noise, shot noise and 
black body radiation) that are fundamentally different into a single expression. It is additive 
so that it is added to the signal, it is white meaning that it has the same power for all fre- 
quencies within the bandwidth of the system and it is Gaussian meaning that its real- and 
imaginary values follow a Gaussian distribution. The ratio between signal- and noise power 
is the signal-to-noise ratio (SNR) defined as 


Elrt)"? | E X E 
E[hn(t)?] ^ NoBT,  NoBT, 


where T; is the symbol time, Tp is the bit time, NoB is the noise power, E, and Ey are bit- 
and symbol energies respectively, r(t) is the received signal and n(t) is the AWGN noise. 

A second noise source (called phase noise) originates from the oscillators downconverting the 
signal from the carrier frequency to the baseband. Instead of producing a perfect sinusoid as 
is assumed in (2.1) the sinusoidal output has a small random phase and amplitude fluctuation 


SNR = (2.2) 


r(t) = A[1 + a(t)] cos (wet + ¢ (t)) (2.3) 


Here A is the amplitude of the wave, a(t) is the amplitude noise, ¢ (t) is the phase noise 
and we is the carrier oscillator central angular frequency. The amplitude noise can usually 
be neglected (because of amplitude limiting mechanisms) while the phase component is more 
difficult to handle. The origins and distribution of oscillator phase noise will not be discussed 
here but the most important thing to note is that phase noise will rotate the signal constel- 
lation (e.g. Figure 2.2) in a random manner and if left untreated will completely ruin the 
reception. A detailed discussion on phase noise can be found for example in reference [22]. 
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2.3.2 Fading 


A radio frequency (RF) signal that is transmitted from one location to another can reach 
the endpoint along multiple paths. This is because the signal can be reflected off surfaces 
and thus in a rich environment there is a multitude of copies of the same signal reaching 
the receiver, each with its own delay, phase and signal strength. These different copies of 
the signal can combine constructively or destructively so that the signal power is attenuated 
and we say that we are in a deep fade. Taking the Fourier transform of the channel impulse 
response the response of the channel to different frequencies can be visualized. If there is 
a notable discrimination between different frequencies within the signal bandwidth B, the 
channel is said to be frequency selective. If the channel is similar within the signal bandwidth 
the channel is flat. The coherence bandwidth Be quantifies the bandwidth for which the 
channel can be considered to be flat. So if the signal bandwidth is smaller than the coherence 
bandwidth we have a flat channel but if the signal bandwidth is lower we have a frequency 
selective channel. If the channel is known, that is if the power delay profile is known the 
coherence bandwidth can be calculated, this is done in [20, ch.3]. 


Since the properties of a wireless channel depend on the surroundings which can be highly 
erratic the channel also changes with time. These changes are either due to the movement 
of the transmitter or receiver or because a change in the environment (e.g. a car moves). 
Coherence time Te quantifies the speed of these changes and is a statistical measure of the 
time period over which the channel impulse response is time invariant. If the symbol time T, 
is smaller than the coherence time, the channel is said to be slow fading meaning that each 
symbol sees the same fading for the whole symbol period. If the symbol time is larger than 
the coherence time the channel exhibits fast fading. [20, ch.3] treats the calculation of the 
coherence time by using the Doppler spectrum. 


A frequency selective channel will distort the signal making the entire system response 
(matched filters at the Tx and Rx as well as the channel) non-Nyquist. This will in turn 
cause Intersymbol Interference (ISI) discussed in the next section. 

These effects are summarized in Figure 2.3. The channel studied in this thesis is characterized 
as exhibiting slow and frequency selective fading. 


2.3.3 Intersymbol Interference 


Intersymbol Interference (IST) is observed when the current received symbol is affected and 
distorted by the surrounding symbols. Looking at Figure 2.2(d) it can be seen that if one 
symbol is shifted in time it will start affecting the surrounding symbols. Since in a multipath 
environment each symbol can arrive at multiple time slots and/or shifted within its own it is 
apparent that ISI can ruin the reception if not properly counteracted. Two other mechanisms 
that can introduce ISI are symbol timing errors and non-Nyquist filter pairs at the Tx and Rx. 


The above effects can be seen by looking at the channel impulse response which can be made 
to include Tx and Rx circuitry, that is filters (p(t)) and analog/digital converters for example. 
An example of a channel impulse response that causes ISI can be seen in Figure 2.4. To see 
how this channel impulse response causes ISI it can be convolved with a symbol sequence sn. 
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Slow and Fast and 


Freq. Freq. 
Bs>>Bc | Selective Selective 
Fading Fading 
Slow and Flat Fast and Flat 
BOSSES Fading Fading 
Ts««Tc Ts>>T¢ 


Figure 2.3: Channel characterization. There are four different kinds of channels based on if 
the signal bandwidth is smaller or larger than the coherence bandwidth and if the symbol 
time is smaller or larger than the coherence time. 


Tn = p» hi8n—i + Mn = hgSn-5 + p» hiSn—i + Nn (2.4) 
E 125 
bu hs 


Figure 2.4: An example of a channel impulse response that would cause ISI. Tap marked by 
hs is the central tap seen in (2.4). The rest of the impulse response will then contribute to 


ISI. 


Here n, is the noise and h; is the channel impulse response. The first term in the equation 
above now represents the desired symbol, delayed and with some gain hs. The second term 
is some undesired contribution from the surrounding signals multiplied with some channel 
gains, this part represents the ISI. 


2.4 Microwaves 


Frequencies higher than 6 GHz are considered to be unusable for mobile communication 
because of increased shadowing by obstruction during propagation [40]. These frequencies are 
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thus more useful when there is a LoS path between transmitter and receiver!. Another issue 
with these high frequencies is their increased attenuation. The received power is proportional 
to the square of the wavelength due to the electrical effective area of the receive antenna, as 
can be seen from the free space propagation equation [20, eq. 2.7] 


VGA 
4rd 


2 
P, 


B 


(2.5) 


where G is the product of the antenna gains, P, is the received power and P; is the transmit- 
ted power, d is the distance between the antennas and A is the wavelength of the carrier. This 
means that going from 0.3 GHz to 30 GHz introduces 10.000 fold attenuation for the same 
setup. This can be compensated for in microwave links using highly directional antennas 
(with large G). 


Another more serious phenomenon that is encountered when going up in frequency is the 
attenuation caused by rain. Figure 2.5 shows the attenuation per km in dB as a function 
of frequency and rain rates for a horizontally polarized signal. The attenuation will differ 
somewhat for different polarizations. For links operating at high frequencies even a short hop 
length can experience outage when there is heavy rain. A common requirement for microwave 
links is that they are operational 99.999% of the time (i.e. 5-9's availability) which corresponds 
to a downtime of less than 5 min per year. The link then needs to be able to cope with all rain 
conditions that appear on average for more than 5 min a year. Heavy rain that statistically 
only is seen for 5 min/year for a given location is commonly referred to as 5 min rain. For a 
certain location this sets an upper limit on the distance over which a microwave link can be 
reliably operated. For Europe the 5 min rain rarely exceeds 50 mm/hr while in south Asia it 
can exceed 100 mm/hr. 


2.5 Dual Polarized Systems 


Electromagnetic waves travel in 3D space along a straight line with orthogonal and oscillating 
electric and magnetic fields (plane wave approximation). The condition is that the fields are 
orthogonal to each other as well as the direction of propagation. A plane that is perpendicular 
to the propagation direction will thus contain all of the electromagnetic fields so that it is 
sufficient to describe them within that plane. For a given wave either of the field components 
(magnetic or electric) can be oscillating for example along the vertical axis requiring the other 
field component to oscillate along the horizontal axis. 

Since most waves of interest for microwave links travel along the earth surface one can speak 
of a vertically- or horizontally polarized waves meaning that the electrical field is oscillating 
in the vertical or horizontal direction respectively. Naturally the wave can be polarized along 
any direction but the important thing is that each polarization can be broken down into 
these two components that form an orthogonal basis. A dipole antenna is a good example of 
a polarized antenna since it will only transmit and receive one polarization in the ideal case. 
Having two ideal polarized antennas rotated so that they are orthogonal to each other results 


l'Today's microwave links operate in frequency bands between 6 and 38 GHz 
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Figure 2.5: Signal attenuation caused by rain for different frequencies for horizontally polar- 
ized signals [12]. 


in two independent channels doubling the capacity of the system for a high SNR. A real sys- 
tem will always experience imperfections leading to crosstalk between polarizations. This is 
due to the fact that antennas can not be perfectly polarized and channel conditions may be in 
such a way that they treat the polarizations differently (e.g. rain) resulting in leakage between 
the channels. This leakage can be quantified using the channel cross-polarization discrimi- 
nation (XPD) factor. It describes how much power from one polarization leaks into another 
polarization thus reducing the systems ability to separate between the two polarizations. It 
is defined for the vertical- and horizontal components respectively as 


E {|v v|} E(|ha,11?) 
E {Ihave} = 


where for example hy; 1 is the flat channel impulse response between the vertically polarized 
Tx and horizontally polarized Rx. The higher the XPD the better isolation between different 
polarizations. Here these two XPDs are combined into one using the assumption that the 
they are the same for both polarizations, i.e. XPD = XPDy = XPDy 


XPDy = (2.6) 
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Chapter 3 


MIMO Systems 


Multiple-input-multiple-output (MIMO) systems are those having (as the name suggests) 
more than one input and output, here this is taken to mean that both Tx and Rx have multiple 
antennas operating over the same bandwidth. There has been a great interest in research into 
wireless MIMO technologies for the last decade which was sparked by ground breaking work 
done by Winters [50] and Foschini [16] in 1992. These papers predicted remarkable spectral 
efficiencies for wireless communication systems with multiple transmitters and receivers. As 
time progressed interest would grow since the value of spectrum seemed to have no limits, this 
was reflected in the radio spectrum auctions in for example Germany and Britain during the 
height of the dot-com boom in 2000. There companies bought relatively small bandwidths 
for 100s of millions of pounds or Euros. The ability to squeeze even more bits through using 
smaller portions of spectrum became very valuable indeed. The interest subsided though 
slightly after the dot-com burst in 2001 but the technology never lost its appeal and now 
there are quite a few standards incorporating MIMO (e.g. LTE, WiMAX and IEEE 802.11n) 
and quite a few products using it. It can be said that MIMO is a proven technology with 
an abundance of published material, a good overview of MIMO can be found for example 
in [40]. A pre-condition for the great majority of this work is that the transmitter and/or 
the receiver are in a rich scattering environment unlike the line of sight links examined in 
this thesis. For a mobile scenario where the handset is surrounded by objects reflecting the 
signals from many different directions without any path dominating, the technology is at its 
best. In fact, the LoS MIMO path between multiple Tx’s and Rx’s have for practical mobile 
applications been considered useless as the difference in path lengths is too short to reach 
orthogonality between the received signals. 

Moving the higher frequencies, e.g. microwaves, and requiring a fixed setup solves this prob- 
lem, enabling MIMO to work even with only the direct paths present, but at the expense 
of larger antenna spacings as well as requirements on the distance between transmitter and 
receiver. It should be noted that the indoor environment can also be suitable since only short 
distances are needed. LoS MIMO is a fairly recent research topic, papers focusing on the 
indoor environment [25, 27] and outdoor environment [18, 14, 39] have been published on 
LoS MIMO over the years but the main attention has always been on the classical MIMO 
multipath environment. Excellent work on LoS MIMO has also been done by Bohagen [5]. 


Besides spatial multiplexing, a second advantage of MIMO systems is that of adaptive di- 
rectivity called beamforming or interference cancellation so that an array of non-directive 
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antennas can be electrically controlled to have a large gain in some direction while cancelling 
others (containing an interfering source for example). Diversity gain can also be obtained 
from the added antennas. These techniques are called smart antenna techniques and are 
discussed in [4, 23, 36]. 


3.1 Narrowband MIMO 


A simple schematic view of a MIMO system can be seen in Figure 3.1. Here data enters the 
MIMO transmitter and is processed and distributed between the Tx antennas. The signal is 
transmitted over a channel described by the channel matrix H. The signals then arrive at 
the Rx antennas combined, it is then the job of the MIMO receiver to separate and process 
the incoming signal streams to reproduce (hopefully) the same data stream as at the Tx. 


MIMO | Data 


Figure 3.1: A block diagram of a simple MIMO system with N transmit antennas and M 
receive antennas. 


Narrowband channels are called narrowband since they treat all the frequency components 
within the bandwidth of the signal the same (symbol duration is much larger than the multi- 
path delay spread T; >> Tm), that is the frequency response is flat. That is to say that the 
received signals are only attenuated, delayed and phase shifted which means that each chan- 
nel gain can be described (neglecting delay) by a single complex number. The whole channel 
which is made up of multiple paths can then be specified by a channel matrix H. Having N 
transmit antennas and M receive antennas and Lg symbols to be transmitted through each 
Tx antenna the system is described by 


Ti occ UM S11 c SiN hi c: him 
fig oe FEQ SLg1 sie Sign hwi © hwM 
nj DM 
+ (3.1) 
NLgl *** NLgM 


where rim is symbol number i from receive antenna m, Sim are the transmitted symbols, Anm 
are complex values describing the channel between transmit antenna n and receive antenna 
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m and finally nim represents the Gaussian noise. Using matrix notation (3.1) can simply 
be written as r = sH +n. The scalar values of the channel matrix H can be written as 
hag = Oamexp(j8ngs) where anm and Onm are the amplitude- and phase of the channel 
respectively. 


The problem of the receiver is to counteract the effects of the channel, one way to do that is 
to find a matrix B so that H.B—IifN«-MorB-.H-—Iif N > M. Each Rx antenna 
receives symbols from all Tx antennas so it is the job of the receiver to separate the data 
streams to obtain an estimate of what was actually transmitted. Assuming receiver channel 
knowledge, that is that H is known at the Rx, B can be found by using the Moore-Penrose 
Matrix Inverse 
Hygp gH ; 
pes d d mol (3.2) 
H'(HH") if N<M 


This can then always be done given that the inverse exists which is always the case if H is full 
rank. The rank of H, denoted Ry, defines how many linearly independent columns a matrix 
has which again corresponds to the number of achievable independent parallel channels that 
can be obtained. The rank always satisfies Ry < min(N, M) meaning that the number of in- 
dependent channels obtainable from H is always equal or less than the lower antenna number 
at the Rx or Tx. The Pseudoinverse in (3.2) will enhance noise! and is not the optimal way if 
there is some prior knowledge about the noise distribution and power (e.g. MMSE estimate °). 


When the channel is known at the receiver or transmitter it is called channel state information 
at the receiver (CSIR) or channel state information at the transmitter (CSIT) respectively. 
There are three basic schemes that are used to obtain CSIR and these are: Training based, 
semi-blind and blind channel estimations. Training based methods [45] send training symbols 
(symbols that are known at the Rx beforehand) and use them to compare what was received 
and what was sent to calculate an estimate of the channel. These training symbols can then be 
sent periodically in certain time slots within the symbol stream to keep track of the changing 
channel. Blind methods (e.g. [52]) are based on the fact that the Rx knows the modulation 
scheme and thus has some knowledge about how a transmitted signal should look like. It 
can thus compare the shape of the received symbols to what could have been transmitted. 
Since what was actually transmitted is unknown these methods can be slow to converge and 
may suffer from some phase ambiguity even after convergence. Semi-blind methods use some 
combination of training based and blind methods. The training data is sent only in the 
beginning [37, 10, 24] or is embedded within the data symbols called superimposed pilots 
(SIP) [8, 9, 48, 47, 26, 35]. They can get rid of the phase ambiguity as well as improving the 
convergence speed compared to blind methods. 

CSIT is then usually obtainable by the Rx sending information about the channel to the Tx 
or the Tx using channel reciprocity where applicable. 


¡Since when removing the channel affects the noise is not considered. 
21f N > M then B = (HE H + 621) ! HP. where 62, is the noise variance. 
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3.1.1 Singular Value Decomposition and Condition Numbers 


As mentioned before not all environments are suitable for MIMO transmission. If for example 
the channel has one dominating path such as in a LoS environment then the capacity of the 
MIMO system should suffer. To see this we need to apply some mathematical method from 
linear algebra, namely singular value decomposition (SVD). For every N x M matrix H one 
can write 


H = UZV" (3.3) 


where U and V are N x N and M x M unitary matrices ? respectively. Y is a M x N 
diagonal matrix of singular values c; of H. If V is now known at the Tx (CSIT) and U 
at the receiver (CSIR) then the symbols at the Tx can be multiplied by UH before being 
transmitted and the received symbols multiplied by V after being received. This parallelizes 
the MIMO channel into Ry independent SISO channels with gains 01, 02, ,... ,0Ry (see 
Appendix B for further discussion). This can be seen in Figure 3.2. Since now these singular 
values represent the gain of each SISO channel the best channel matrix is the one that has 
equal gains or in other words the singular values should be identical. The larger the ratio 
between these singular values the worse the MIMO channel is. To describe the theoretical 
performance of a MIMO channel one needs to know all the singular values. It is though clear 
that it is always the smallest singular value that will affect the BER the most. In this thesis 
the condition number is introduced as a measure on how good a MIMO channel is. The 
condition number is the ratio between the largest- and the smallest singular value. When the 
channel is perfectly conditioned the condition number is one but it will grow as the channel 
condition worsen. The condition number is defined as 

je (3.4) 

min(0;) 

The number has no mathematically rigorous coupling to BER or penalty but a rule of thumb 
is that if the condition number is large the penalty should be close to the square of the con- 
dition number. That is the SNR needs to be raised by x? to achieve the same performance as 
if the condition number was one. One has also to remember that penalty is usually expressed 
in dB while condition number is in a linear scale. 


As an example lets consider two 2x2 MIMO channels H; and Hs given as 


1  em/2 1 1 
m - [t i | and H= | i " (3.5) 


Doing the SVD so that Hy, = UX VË and H5 = U3X,VH we get the following singular 
value matrices 


m= [Vo] ua me[?] 


The singular values for the first channel are 011 = 012 = V2 while the singular values for the 
second channel are 621 = 2 and 099 = 0. The condition numbers are thus 


3A unitary matrix H satisfies H"H = HH"! =I 
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Figure 3.2: Using SVD the MIMO channel can be parallelized into Ry independent SISO 
subchannels with gains ox. 


K1 = v2 
y? (3.6) 


K2 = 7 œ 


The second channel is rank deficient (Ry = 1) and does in fact not support MIMO transmis- 
sion while the first one is full rank (Ry = 2) and is ideal for MIMO. This can also easily be 
seen if we try to invert H4 and Ha since Ho is singular and not possible to invert. 


3.1.2 Channel Capacity 
The capacity of a SISO system is given by the Shannon-Hartley theorem 


C = log, (1 + y) (3.7) 
where y is the average SNR. The capacity of a MIMO system was derived by Telatar in 1995 
[44]. The result is 


E min(N,M) - 
C = log; [det (Iu j 7 HH")| = M dw (1 J Ta) (3.8) 


where A; is a eigenvalue of HH". The eigenvalues and singular values are related by A; = a. 


Looking back at the matrices defined in (3.5) the capacity can thus be calculated using the 
singular values obtained in the previous example for a given SNR. Setting the SNR to 20 dB 
the capacity becomes 


2 
Cr = Y og, (1 +100 - v2) ~ 15.3 [bps/Hz] 
i—1 
Cz = logy (1+ 100.2?) + log, (1+ 100-0?) ~ 8.6 [bps/Hz] (3.9) 
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As can be seen the first channel has almost twice as high capacity compared to the second 
one. This is because H> has rank 1 and only supports one data stream. If the SNR has a 
lower value the capacity gain of case 1 in the above equation will be reduced. Thus MIMO 
works best at high SNR which is most often readily available when there is LoS. 

It is also noteworthy to mention that if the Tx knows the channel (CSIT) he can allocate 
more power to the better channel thus increasing the capacity (e.g. waterfilling). Assuming 
that the Tx in case 2 allocates double the power to the better channel and no power to the 
other one and the SNR has the relatively low value of 4 dB the result becomes 


Ci 


2 
Y log, (1 125. v2) = 5.2 [bps/Hz] 
i—1 
Cy = log) (1+2.5-2-27) &: 44. [bps/Hz] (3.10) 


3.2 Line-of-Sight MIMO 


A properly designed MIMO system in a rich multipath environment will have a low correla- 
tion of phase and amplitude between different receive antennas. Replacing the rich multipath 
environment with one dominating LoS path will correlate the channels so that there will al- 
ways be the same phase offset and amplitude between the Tx- and Rx antennas pairs over 
time. If the antennas are placed closely together the channel matrix will become a matrix 
with all elements identical which is singular reducing the capacity to a SISO system (3.7). 


In LoS the phase of each channel gain hnm is a function of the geometry and carrier fre- 
quency but is not governed by a statistical distribution as in non-LoS MIMO. By varying the 
separation between the transmitting and receiving antennas a static phase shift between the 
antennas can be introduced thus effecting the channel matrix. This is illustrated in Figure 
3.3 where the different paths are marked ham. The phase of the received signal stream will 
depend on the length of the path (and carrier frequency) so different paths can have different 
phases. The normalized channel matrix can then be expressed as 


hi ©- him exp (jkd11) ++» exp(jkdim) 
H= HOS : = : ET : (3.11) 
hyi ©- hww exp(jkdwi) ~ exp(jkdym) 


where k = 27/X and dnm represents the physical distances between the antennas. 


3.2.1 Optimal Antenna Separation 


As mentioned earlier the LoS channel is described by the geometry (in units of wavelengths) 
of the antenna placements. Finding the optimal configuration is thus a geometrical problem. 
To solve it some preferred antenna arrangement must be assumed. Here the antennas are 
considered to be organized in linear arrays where the Tx- and Rx antennas are separated by d; 
and d, respectively as shown in Figure 3.4. The total length of the arrays are then (N — 1)d; 
and (M — 1)d,. The Tx and Rx antennas can be tilted in the plane of transmission by 0, and 
0,. respectively and the Rx array can be rotated out of plane by $, as is shown in the right 
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Figure 3.3: A schematic view of a 2x2 LoS MIMO system. There are four channels hij 
describing the phase and amplitude of each path. 


side of Figure 3.4. Considering the origin to be at Txo and denoting n and m to be the index 
of Tx and Rx antennas respectively, starting with zeros and going up to N — 1 and M — 1, 
then doing the derivation similar as in [5] the vectors describing the antenna placements are 
given by 


at = nd, sin(6;)n, + nd; cos(0,)n; 
a”, = [D + md, sin(6,) cos($,)]n, (3.12) 
+md, cos(0,)n; + md, sin(0,.) sin($,.)ny 


where n;, n, and n; are unit vectors along the x-, y- and z-axis respectively. To calculate the 
distance between antenna m and n the absolute value of the difference of the vectors above 
can be taken. This corresponds then to the Euclidean distance between the antennas 


dnm = la”, —a! | = [(D + md, sin(0,) cos(¢,) — ndisin(0;))? 
+ (md; sin(6,) sin(¢¢))? + (md, cos(8,) — nd; cos(6,))? |? (3.13) 


Using the approximation v A? + B = A+ A if A >> B, which is true for the equation above 
since the distance D is much larger than the antenna spacing, then (3.13) can be written as 


dnm œ~ D + mdrsin, cos $, — nd; sin 0, 


1 
+ 3D (md, cos 0, — nd, cos &y + (md, sin 6, sin br)? 
= D + md, sinb, cos ¢, — nd; sin 0; 
1 
+ 3D [m?d? cos? 6, + n?d? cos? 0, — 2mndid, cos 0, cos 0,. + m?d? sin? 6, sin? or] 


Here another simplification has been made which is D + md, sin(0,) — nd; sin(0,) ~ D for the 
denominator which is true since the distance between the Tx and Rx is much larger than the 
antenna separation and is also true for small angles. 
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Figure 3.4: Tx antennas and Rx antennas organized in an array. The separation between the 
Tx and Rx antennas is d; and d, respectively, while 0, and 6, represents the angles of the 
antenna arrays with respect to the z axis. D is the distance between the antenna arrays. 


Now that the distances between the antennas have been calculated a criteria for the optimal 
antenna separation should be found. One way to accomplish this is to note that the optimal 
separation should give the highest capacity or in other words should maximize (3.8). This 
corresponds to maximizing the product of the eigenvalues A;. It is known from linear algebra 
that trace{W} = >> \; = NM for a normalized unit power matrix H if 


v= HH!, N<M 
— | HEH, N>M 


Then it is easy to show that (3.8) is maximized by A1 = Ag = + +- = constant. This is obtained 
if H has orthogonal rows for N < M or orthogonal columns for N > M. Defining the rows 
of Hzos as hn, orthogonality between them can be expressed as < hy, hj »,7;— 0. This 
condition can then be used to produce the expression for the optimal antenna separation (if 
N <M) 


M-1 
< Tij; hi Dngi = > exp (Jk(dum dim)) 


M-1 
d, cos(0,. 01),. 
= b» exp (20% coste a ae nym) 


m=0 


exp E — n)d; sin(0,) + 5p — n)?d; cos*(6)]) =0 


4 sin (kd;d, cos(0,) cos(0+) (i — n)M/2D) =i 
sin (kd;d, cos(0,.) cos(0,) (à — n) /2D) 
AD 
MP dip M cos(0,) CE ety) 


where K is a positive odd number usually chosen to be one since that gives the smallest 
optimal antenna separation. Doing a similar derivation for the case N > M would yield 
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the same expression only with M replaced with N. Defining V = min(M, N) a general 
requirement on the antenna separation is given as 


AD 
B K 
V cos(0,) cos(0,.) 


Looking at the optimal antenna spacings given in (3.15) it can be seen that the separation 
increases as the distance D increases and that it decreases as the carrier frequency fe increases 
(since À = c/ fe). From an installation point of view it is often practical to chose 0; = 0, = 0 
and K — 1. This is done in Figure 3.5 where the optimum antenna separation is plotted 
against distance and for different carrier frequencies. 
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Figure 3.5: Optimal antenna separation for a 2x2 SS LoS MIMO system, plotted against 
distance for a few carrier frequencies. If the separation at both the Tx and Rx is chosen to 
be the same d; — d, — d then the values on the y-axis correspond directly to d. 


Other geometries of the antenna setup at Tx and Rx can be chosen as well. One such ap- 
proach is to arrange the antennas in an equilateral rectangle (minimum of 4 antennas needed 
at Tx and Rx) also called uniform planar array. The calculations yield almost the same result 
as for the linear antenna array as is shown in reference [6]. 


There might be a few different reasons why the antennas are not placed with the optimal 
separation. The most common ones are perhaps that the required separation might be too 
large to accommodate at the site and that all transmission systems have some bandwidth, it 
is only possible to optimize for one frequency. Defining the parameter y as 


n= did, 
(didr opt 


to quantify the deviation from optimality. Choosing d; = d, = d, then (3.16) reduces to 


(3.16) 
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n = d/ dopt- The condition number as a function of ņ for a 2x2 LoS MIMO is plotted in 


Figure 3.6. The condition number is at its lowest at y = VK where K = 1,3,5,..., the 
condition number approaches infinity for K = 0,2,4,.... 


Condition Number Vs n 


Condition Number 
a 


Figure 3.6: The condition number as a function of 7 defined in (3.16) for a 2x2 LoS MIMO 
system. 


The condition number for a 2x2 LoS MIMO is derived analytically in Appendix A.1 and is 
given by 


2 + [2 + 2cos(12)] 1? 


2 — [2 + 2cos(12)] 2 Ven 


K = 


3.2.2 Vector Visualization for LoS MIMO 


An alternative way of understanding a LoS MIMO system would be to use vector notation. 
For a 2x2 LoS MIMO system it can be seen that for the optimum case there should be a 
90° phase shift at the carrier between Tx antenna n and Rx antennas m ¥ n (because of the 
condition < h;, hn >iżn= 0). This can be visualized graphically using a vector representation 
of the signals. Consider first Figure 3.7a) where the Tx- and Rx antennas are phase locked 
and there is a 90° phase shift at the carrier frequency between paths dnn and dnm where 
n z m. Both receive antennas receive a signal that is a combination of both transmitted 
signals only with different phases. By rotating the received signals and combining them the 
combined signal streams can be separated. This is true even though the transmitters are not 
phase locked (Figure 3.7b-c). If however the phase difference between the direct- and diagonal 
paths is less than 90° as it is in Figure 3.7d the resulting signals will have less power and are 
thus more susceptible to noise. This corresponds to a suboptimal situation which can also be 
seen from the high condition number of the LoS channel matrix. 
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Figure 3.7: Separating the signals streams in a 2x2 LoS MIMO. All figures assume that the 
distances d;; are a multiple of the wavelengths so that the received signal has the same phase 
rotation as the transmitted one. 90° difference between the direct- and diagonal paths for 
a-c, 45? for d. a) The Tx signals are phase locked (0? phase difference between them). b) 45? 
phase difference between Tx signals. c) 90? phase difference. It is interesting to note that 
there is no signal strength at Rx antenna 2. d) 0°. 


Figure 3.7c) is an interesting special case. There Rx» receives no power while Rx; receives 
double, this happens because the signals combine constructively at Rx; while they combine 
destructively at Rx2. Having a constant modulus signal constellation such as 4QAM the 
symbols will always have the same power but will instead be rotated differently. If Tx» 
transmits a symbol with a 90° phase shift compared to Tx, as is depicted in Figure 3.7c) Rx? 
receives no power. But if the receivers follows the simple procedure of rotating the received 
symbol back 90? and summing up between the branches the correct symbols will still be 
recovered. 


3.2.3 Grating Lobes 


A third approach for understanding LoS MIMO would be through grating lobes. Having 
two spatially separated transmitting antennas operating in the same bandwidth and carrier 
frequency will create gain patterns in the far field called grating lobes. For the setup in 
Figure 3.8 the receiver will see different signal strengths depending on the physical angle @ 
and electrical angle a for a fixed Tx antenna separation and frequency. We will start with 
the case when both Tx’s are transmitting the same signal which will result in a static gain 
pattern. 

Calculating the distances dı and d» using the law of cosines one gets 


= d?/4+ D? — d D cos(6) 
di = d? /4 + D? + dD cos(6) 


The gain (relative to a single Rx antenna) can also be expressed as 


(3.18) 
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1 Tse 
G(6,a) = z- : (sin(wt + dık) + sin(wt + dak + a))? dt 
= cos(dik — dək — o) + 1 


= cos(kAd+a) +1 (3.19) 


where Ad = də — di. 


di 


Figure 3.8: Two Tx antennas both transmitting the same signal. The observed signal strength 
or gain will depend on the angle between the Tx and Rx. 


The grating lobes should have the same shape independent of D if the receiver is situated in 
the far field. For this to be true Ad needs to be constant if @ and d; are constant. 


(da —d1)? = (di-dj-2did;) 


2 4 
= E + 2D? — E + 4D* + 2d?.D? — 4d? D? cos?(¢) 
ee d; 12D? 2p? dt /4 + 2d? D? — 4d? D? cos? (¢) 
2 2D? 
2 2 di 
= d¿cos*(p) — 16D? (3.20) 


where the approximation VA? + Bx A+ E if A? >> B was used in the third step. For 
D >> d (3.20) can be written as 


Ad ~ d;cos(¢) (3.21) 


For two transmit antennas separated by a certain distance, (3.19) describes the power gain 
of the received signal for an observer in the far field. Envisioning an observer moving around 
the antenna array in the far field and measuring the signal strength it can be seen that as the 
physical separation in wavelengths between the transmit antennas grows the faster the signal 
strength will vary. Combining (3.19) and (3.21) it is easy to see how this is true 


2 
G(¢, a) = cos(kAd + a) + 1 = cos (Fa cos(@) + a) +1 
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As the ratio d;/A increases then the gain expressed above will vary faster as $ is varied. This 
can also be visualized looking at Figure 3.8, the larger d; is the larger the change in Ad for 
some Af. If now the transmitter would like to put the receiver in a null gain it can do so 
by varying a. a represents the electrical angle that is the phase of the signal and thus a full 
180? degree phase offset is required to go from full gain to a null. The effects of changing « 
and varying the frequency (or antenna separation) can be seen in Figure 3.9. 
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Figure 3.9: Gain pattern (grating lobes) for two antennas separated by d;. The gain varies 
faster with angle as the antenna separation increases in wavelengths. 1) a = 0 and d; = 2X. 
2) a = 0 and d; = 8X. 3) a = 1 and d; = 24. 4) a = s and d; = 8X. 


As the derivation in section 3.2.1 showed, the separation between Rx antennas and Tx an- 
tennas (denoted d, and d, respectively) should be in such a way that there is a 90? phase 
shifts between the direct- and diagonal paths. Setting the antenna separations so that (3.14) 
is satisfied and plotting the gain each Rx antenna sees while rotating the Tx antenna array 
an interesting pattern emerges. This can bee seen in Figure 3.10. When Rx» has a deep fade 
Rx; has maximum gain and vice versa. In LoS MIMO we are usually interested in spatially 
multiplexed systems so that each Tx antenna is transmitting different symbol streams. If 
the symbols on each Tx antenna have the same amplitude (by coincidence or because the 
modulation is constant amplitude) then their phase difference correspond to o. Since o only 
shifts the pattern seen in Figure 3.10 the situation can arise (for a 2x2 SS LoS MIMO) that 
one antenna receives no signal while the other receives a combination from both Tx antennas. 
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That situation is depicted also in Figure 3.7c and it is interesting to note that the signal 
streams are still easily separable at the receiver. 
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Figure 3.10: Gain pattern for an antenna arrangement satisfying (3.15). When one Rx has a 
null gain the second Rx will experience a peak gain. 


3.3 Wideband MIMO Model 


The discussion so far was limited to narrowband MIMO or flat channel MIMO, that is channels 
that could be described by a single complex number. Thus it was easy to express the MIMO 
channels in matrix form H so that the received signal was r = sH +n. That can also be 
expressed in a series form as 


N 
n=1 


Here m represent the index of the Rx antenna and j is the discrete symbol timing. Now a 
wideband channel will act as a filter so that the gains in (3.22) need to be replaced by filters. 
This can be accomplished by 


N It 
rm[j] = `y y hnmli, jlsnlj — i] € 


n=1 i=— Li 


N 
rm = > hnm * Sn 
n=1 


where hy represents the channel filter between Tx antenna n and receive antenna m of 
length 2L; + 1 and s, is the vector representation of the symbols sent from Tx antenna n. 

Expressing frequency selective MIMO channels in matrix form is a little bit more difficult 
and it requires some modifications to the channel matrix. In the most general case a N x M 
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MIMO system has NM channels where each one of them is made up of a certain number 
of multipaths. The full channel is then considered to be made up of NM filters that are 
uncorrelated in the most general case. Equalizing such a channel is a much harder task than 
equalizing SISO channels since NM independent filters will be required in stead of just one. 
Transmitting Lg symbols (also called block length) the channel matrix can be expressed using 


Toeplitz matrices 4 such as 
hob dix ee y — a OT 
"n hnm[1] hnm [0] ps Bus su 0 


where the Lp x Lg matrix W,,,, describes the channel from transmit antenna n to receive 
antenna m and 2L, + 1 is the length of the channel filter. Expressing the full MIMO channel 
in this way gives 


fies nhi dl 
w- ro: 
EN eni Pyu | 
Defining the signal matrix s as 
s = [s1,...,Sy] 
Sn = lsali suba 


where snpli] is the ith symbol transmitted from Tx antenna n. Here s will be a 1 x NLg 
vector. This enables the noise free received signal to be expressed as 


r= sọ (3.23) 


where r is the received signal and will have the dimensions 1 x MLpg. The matrix Vn», 
arranges the channel filter taps in a convolutional matrix form. For a SISO channel there is 
only one filter so YW ; multiplied by the signal vector would represent the filtering operation 
of the channel. For a MIMO system there are multiple channel filters that operate on the Tx 
signal streams which are added up at the receiver. This may be conveniently expressed in 
matrix form by (3.23). 


3.3.1 Channel Capacity 


The channel capacity was defined for a narrowband MIMO channel in Section 3.1.2. For a 
frequency selective channel the capacity becomes frequency dependent. For a SISO channel 
the capacity becomes [33] 


C= P logy (1+ HD) df 


4A Toeplitz matrix is a square matrix with equal values on its left-to-right diagonals. 
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where H(f) is the transfer characteristic of the channel (including Tx and Rx filters). A 
MIMO channel is described by a channel matrix so that the capacity can be expressed as [34] 
min(N,M) A, 
c= J m (in 
i=l nb 


where A;(f) is the ith eigenvalue of the positive definite matrix C(f) = HE(f)H(f). Here 
H( f) represents the frequency dependent narrowband channel matrix. 
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3.4 The LoS MIMO Channel 


The LoS-MIMO channel differs from the non-LoS MIMO channel by the fact that there is one 
dominating path so that the incoming phases are dominated by the geometry of the channel 
and not scattering. Besides the LoS-path there can be secondary paths caused either by 
atmospheric scintillation or a reflection off the ground. This results in a frequency selective 
channel. A common model to use is the so called two ray model where the channel is described 
using one direct path and one path that reflects off the ground and is received with a random 
phase ¢ and a relative amplitude of 0 € b < 1. The two paths have some delay difference 7 
that is commonly set to 6.3 ns as in Rummler's model [3]. The impulse response of such a 
channel is 


A(t) = S(t) +bå(t-— r)e? S 
H(s) = £(h(0)1-be 7939 (3.24) 


Remembering that s — jw it can be seen from (3.24) that the channel introduces periodic 
notches in the frequency spectrum which are illustrated in Figure 3.11. The channel spans 
some bandwidth and will thus have a frequency response depending on where it is situated 
in the spectrum of Figure 3.11. If one of the notches happen to be within the bandwidth the 
channel will suffer from severe fading and ISI. The depth of the notches is given (in dB) by 
B, = 20 log; 9(1 — b). 
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Figure 3.11: Periodic notches for the Two Ray Model. 7 = 6.3 ns, ¢ = 0 and b = 0.9 in 
(3.24), 


A LoS 2x2 SS-MIMO system with four frequency selective channels, denoted hy, is illustrated 
in Figure 3.12. By varying the phase $44 of the ground reflection the properties of the 
channel filters hnm are changed. Each channel can act as a lowpass- highpass or notch filter 
for example depending on ¢nm. 

The N x M LoS MIMO channel is made up of NM such filters, each with its own realization 
of Ónm. 
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Figure 3.12: A LoS 2x2 SS-MIMO system with four frequency selective channels. 


3.4.1 Mast Swing 


A common installation arrangement for backhaul links is that they are mounted on a tall 
mast along with the mobile base stations which are placed on the top with the backhaul 
antenna disks a little bit below. In heavy wind the mast might rock back and forth resulting 
in periodic movement of the antennas. As the antennas move, the channel between Tx and 
Rx will change since the channels are partly described by geometry. As long as the Rx does 
not move out of the main beam of the Tx antenna the signal strength between any pair of Tx- 
and Rx antennas should be relatively unaffected. By moving the mast the channel notches 
can be seen across the bandwidth of each channel forcing the STE to adapt. Another effect 
is that the movement at such high carrier frequencies will cause some differential frequency 
offset at the Rx. Since all real systems have to have some carrier recovery this is usually not 
such a large problem. 


3.5 Dual Polarized MIMO 


Dual polarized MIMO systems make use of two orthogonal polarization states, often described 
as vertical and horizontal polarization. The interference between the two states is called cross- 
polarization discrimination (XPD) (defined in (2.6)) and depends both on the transmission 
channels and the polarization mismatch between Tx and Rx. XPD is usually expressed in dB 
and the larger it is the more isolation there is between the states. 


The channel can lower the XPD and there have been quite a few papers on how the channel 
treats the polarizations, e.g. [2, 30]. The results indicate that even in a scattering environ- 
ment the XPD should be high (at least 5-15 dB) and for a non-scattering LoS channel it 
is even higher. Thus DP MIMO schemes can provide high rank MIMO channels with good 
condition numbers. References [15, 42] show the benefits of using DP MIMO in a Ricean 
fading environment when the antennas are not configured for LoS. This is though not the 
optimal comparisons for the purpose of this discussion since the antenna placements for a 
SS MIMO system should be optimized for the LoS environment. In Ref. [7] this is done 
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with the result that SS-SP MIMO outperforms its DP counterpart in low SNR because of the 
added 3 dB array gain. The difference is reduced when the SNR is high which is common 
in LoS environments. The SS-SP MIMO was also found to be more sensitive to polarization 
mismatch than DP MIMO. Reference [12] shows that the XPD is not lowered significantly 
even for heavy rain which suggests that SS-SP MIMO should have superior performance even 
while it rains because of reduced SNR. A benefit of the DP MIMO system is that the antennas 
may be combined, thus there is no need for spatial separation at the Tx and Rx site. 

The largest benefit for DP MIMO systems is though the relaxed requirements on antenna 
separation, DP antennas can be situated as close together as the installation allows. 


It is interesting to analyze 4x4 SS-DP MIMO systems since these can produce 4 fold capacity 
compared to a SISO system (at a high SNR) while relaxing somewhat the requirements for 
the separation. The narrowband channel matrix can be written as 


hivav hivim hiv2av hiv 2H 
H hinav finan himav him 2H SE E] 
hoviv hav1m hevev hev 2H Ho; Hə 
hogiv həman hnv hann 
yT=aekda ^ aedi yI adii qed dia 
7 ac) *du VI — aes hau Jae kh VI — ae d 
= VA — aei*dn yaeikda VI — aedkdaz Vae) kd22 
ac) *dzi VI — aoei*da Jae k422 VI — aedhdee 


or using the Kronecker product Y as 


ei*dii ejkdia vl-o vya 
H = Hros 9 Wxpp = | ejkd21 ejkd22 e Ja Vl-a 


where a = E hav?) =E hiv?) Here hv mu denotes the gain between the vertically 
polarized Tx antenna n and the horizontally polarized Rx antenna m. o measures the ratio of 
the power for one polarization that is transferred to the other polarization and is constrained 
by 0 « a € 1. In terms of the XPD defined in (2.6) a can be expressed as 


(3.25) 


1 
T XPDSd 
The condition number of (3.25) is plotted against XPD in Figure 3.13. As can be seen the 
condition number drops off quickly with increased XPD. If the antenna placement is subopti- 
mal the system can be seen to be more sensitive to the XPD. Comparing this result to Figure 
3.6 it can be seen that the condition numbers converge to the same value. 


To calculate the condition number of (3.25) analytically, the eigenvalues of HH! should be 
calculated first. The result is (given that 0, = 0, = 0) 
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Condition Number Vs XPD 
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Figure 3.13: The condition number « as a function of XPD for a 4x4 SS-DP MIMO as is 
described in (3.25). This is done for 7 = 0.5, 0.7, 1 (see (3.16) for the definition of 7). The 
dots represent values calculated using (3.27). 


cq = (2- fala)) [1 — cos (xn? /2)] 
ca = (2— fala)) [cos(rn?/2) + 1] 
c3 = (2+ fala)) [cos(ry*/2) + 1] 
ca = (2+ fala) [1 — cos(rn?/2)] (3.26) 


where 7 is defined by (3.16) and fala) = 4Va— a?. The condition number is thus 


k= aa C2, 3, C4) (3.27) 
min(c;, C2, C3, C4) 

Looking at (3.26) and (3.27) there are a few things to note. When a = 1/2 (XPD= 0 dB) then 

C, = co = 0 meaning that the channel can only support 2 data streams. This can perhaps be 

seen intuitively since then there is no discrimination between the polarization states, that is 

the 4x4 SS-DP MIMO has been reduced to a 2x2 equivalent. 


3.6 Penalty 


A perfect system setup, that is a full rank MIMO channel and a system that only introduces 
white Gaussian noise, will have a certain BER for a given SNR. If a target BER is set then 
the SNR that achieves that can be computed. If any imperfections are then introduced in 
the system a higher SNR should be needed to get the same target BER. The difference in 
SNR values is here termed as penalty. It is the penalty in power that has to be paid if this 
imperfection is introduced. For example if there is a penalty of 3 dB present in a system the 


Al 


Tx power will have to be raised by 3 dB to achieve the same BER. 


Assuming Gray coding and a square M-QAM constellation the BER can be calculated for a 


SISO system as 
A 1 3k Es 
n=(1-5)2 y Ex) 


where k = log2(M). For a MIMO system such as is depicted in Figure 3.2 there are M Rx 
antennas ? and equally many parallelized SISO channels with singular values 01, ...,0R;- 
The BER for a MIMO system using SVD is thus given as (similar as in [49, eq. 3.50]) 


i E 1 | 3k E 
P,= Lm 3 (: — Ai) Q ( ime] (3.28) 
The singular values (or gains) c; are calculated for 2x2 SS-SP LoS MIMO in (A.1) and for 4x4 
SS-DP LoS MIMO in (3.26) (since the singular values are the square root of the eigenvalues of 
HH). Selecting the BER P, = 107? and M = 16 (16Q AM) the penalty can be numerically 
calculated for a 2x2 SS-SP LoS MIMO and a 4x4 SS-DP LoS MIMO against 7 and XPD. 
This has been done in Figure 3.14. 


°M represents the constellation dimension while M represents the number of Rx antennas. 


42 


2x2 SS-SP LoS MIMO, 16QAM 4x4 SS-DP LoS MIMO, 16QAM 
20 T T 


Penalty [dB] 
Penalty [dB] 


10 15 20 25 
XPD [dB] 


Figure 3.14: The penalty plotted against 7 and XPD for 2x2 SS-MIMO and 4x4 SSDP-MIMO 
respectively assuming the SVD method. The target BER is 107°. 
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Chapter 4 


MIMO Equalization 


This chapter introduces MIMO equalizers, one of the most important components in MIMO 
receivers. MIMO equalizers, on one hand, work as a typical equalizer to remove the Inter- 
Symbol-Interference (ISI) which is a common phenomenon mainly caused by a frequency 
selective channel, such typical equalization techniques, especially for SISO equalizers, have 
been extensively investigated in the literature [53], [20], [43]. On the other hand, in addition 
to removing the ISI, MIMO equalizers cancel the Inter-Stream-Interference. Figure 4.1 shows 
a schematic view of a 2x2 LoS MIMO system ho; and h12 interferes with hy, and hao» at Rx, 
and Rxg respectively. 


h11 
Tx Rx, 
Txo Rx 
d) h22 (> 


Figure 4.1: A schematic view of a 2x2 LoS MIMO system. There are four channels hij 
describing the phase and amplitude of each path. 


Figure 4.2 shows how a MIMO Space Time Equalizer (STE) uses Constant Modulus Algo- 
rithm (CMA) to remove both the ISI and Inter-Stream-Interference for a 2x2 SS-SP MIMO 
system using 16QAM modulation. As shown in the figure, the received signals are wildly 
distributed before equalization. The right figures show how the transmitted signals are recov- 
ered using equalization. It is noteworthy that the recovered signals are shown after the phase 
recovery loop is turned on otherwise these recovered signals might have a random rotation. 
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Source signals for Antenna2 Received signals without equalization for Antenna 2 Equalized Signals for Antenna 2 
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Figure 4.2: 16- QAM constellation. Left figures: The constellation of the transmitted symbols. 
Middle figures: The received symbols distorted by ISI and Inter-Stream-Interference. Right 
figure: The equalized symbols. 


4.1 Equalizer types 


Equalization techniques fall into two broad categories: linear- and non-linear equalizers. Com- 
pared to non-linear equalizers, linear equalizers have larger noise enhancement which will be 
discussed in the following section but they are more easily understood and simpler to imple- 
ment. Due to the complexity of non-linear equalizers from a practical standpoint, only linear 
equalizers are discussed and compared in this report. We look into MIMO linear equaliz- 
ers consisting of trained algorithms such as Least Mean Square (LMS), and Recursive Least 
Square (RLS) and blind algorithms such as Constant Modulus Algorithm (CMA), Modified 
Constant Modulus Algorithm (MCMA), Simplified Constant Modulus Algorithm (SCMA) 
and their variants. Figure 4.3 shows the categories of equalization techniques. All the equal- 
izers in this thesis use a symbol-by-symbol transversal filter structure which is shown in Figure 
4.4 where wo stands for the central tap. 


4.2 Linear Equalizers 


4.2.1 Linear equalizers structure and noise enhancement 


Figure 4.5 illustrates a block diagram of a linear equalizer. As discussed above, linear equal- 
izers typically suffer from noise enhancement hence they should balance ISI mitigation and 
noise amplification. Since the overall system is seen as a linear system in which the received 
signal can be expressed as an addition of distorted signals caused by ISI and additive white 


45 


Equalizer 


Linear 


i Non-Linear 
Equalizer Equalizer 
Trained/Adaptive Blind 
Algorithms Algorithms DFE MLSE 
LMS | 
RLS CMA 
MCMA 
SGMA Trained/Adaptive Blind 
Algorithms Algorithms 
LMS CMA 
RLS MCMA 
SCMA 


Figure 4.3: Equalizer types, algorithms. 


Gaussian noise (AWGN). This is expressed as 


Lc 
r(k) = Y h()s(k— i) n(k) (4.1) 
iL. 
where r(k) is the kth sample of the discrete received signals at the receiver front end, s(k) is the 
transmitted signal and n(k) is the Gaussian noise. Using the discrete time Fourier transform 
(DTFT) and passing the received signals through the equalizer, the following equation is 
derived 


Y(w) = RW) Heq(w) = S(w) A (w)Heq(w) + N(w) Heg(w) (4.2) 
where the capital letters of the signals and channel response denote their discrete time Fourier 


transform, ISI introduced by channel H(w) can be removed by introducing an equalizer 
Heq(w), in the frequency domain, defines as 
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Figure 4.4: Transversal Filter Structure. 


Noise 


Data Source Channel Equalizer 
H(w) HQ) 


Hew) = az (4.3) 
Yw) = Síw += (4.4) 


Here the last term is the colored Gaussian noise (frequency dependent), and all the ISI has 
been successfully removed from the received signal. The inverse equalizer which has a fre- 
quency response of Heg(w) is often approximated by a finite-impulse or transversal filter as 
illustrated in Figure 4.4 impulse responses are the filter taps w; and the length of the filter 
taps L is expressed as L = 2L, + 1. As can be seen in the figure, the structure looks like a 
convolution operation which convolves the equalizer taps with the received signals in order to 
invert the CIR (Channel Impulse Response) and cancel the ISI. However, if H(w) has deep 
attenuation at any frequency within the bandwidth of transmitted signals, Heg(w) will have 
a high gain in the same range and the noise power range will be greatly enhanced. In this 
case, even though the ISI is removed by the equalizer the system may still perform poorly 
due to a reduced SNR. Thus the goal of equalization is not only to remove the ISI but also to 
prevent enhancing the noise power too much. Conventional linear equalizers process equaliza- 
tion approximately by inverting the channel response and the noise power will be enhanced. 
Equalizers based on the minimum mean square error (MMSE) algorithms have the benefit 
that they are balancing ISI removal and noise enhancement as they are trying to minimize the 
error between the estimated- and desired symbols. Both LMS and RLS are kind of MMSE 
equalizers. 


AT 


4.2.2 MIMO Space Time Equalizer (STE) 


In order to cancel both the ISI and the Inter-Stream-Interference, a MIMO STE has been 
considered. A STE has to be able to not only initialize filter taps but also to update the taps 
iteratively as the channel change over time using some adaptive algorithm. A STE for both 
2x2 SS-SP and 4x4 SS-DP MIMO systems are considered they can easily be extended into 
a MIMO system with N transmitting and M receiving antennas. The channel between each 
Tx and Rx antenna pair is considered to be an independent frequency selective channel with 
additive white Gaussian noise (AWGN). 
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Figure 4.6: MIMO Space Time Equalizer diagram. 


Figure 4.6 shows a block diagram of a MIMO STE where the overall system can be modeled 
as a linear system shown in Figure 4.7. Here N and M denote the number of transmitting 
and receiving antennas respectively. h and w denote the channel filters and the equalizer 
filters respectively. The following equation shows the discrete output of the mth receiving 
antenna 


M 
rm(k) = X. s. (k) * hi (k) + nm(k) (4.5) 

m=1 
here rm(k) is the discrete output of the mth receiving antenna at sampling time k and nm(k) 
is the AWGN added at mth receiving antenna at sampling time k. s,(k) is the discrete 
input to the nth transmitting antenna, hn,m(k) is the discrete channel response from the 
nth transmitting antenna to mth receiving antenna. The above equation can be rewritten in 

matrix form as follows 


r=sW+n (4.6) 


where r is the received signal with the dimension 1x M Lp. s and Y are defined in (3.23) and 
(3.3) respectively. The equalized signal can be described by the following matrix 
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y(k) = rW (k) (4.7) 


here y(k) is a vector with M equalized symbols at sample time k, W (k) is the equalizer tap 
matrix that is defined in (4.14), rẹ is the kth row vector of receiving matrix R that defined 
in (4.18). Figure 4.8 shows how the real and imaginary part of the magnitudes of STE fil- 
ter taps looks like after the equalizer converges. The first 43 taps are the coefficients of the 
equalizer wj; and the remaining 43 taps belong to the canceler w21 that cancels Inter-Stream- 
Interference from Tx». 


MIMO Channel MIMO STE 


Figure 4.7: Structure of the MIMO space time equalizer. 


4.3 Trained Equalization 


There are two different equalizer setting approaches, preset and adaptive. For the preset ap- 
proach, the information of the channel response is obtained from the channel measurement or 
channel estimation to set the equalizer filter taps. For the latter approach, instead of estimat- 
ing the channel, equalizer taps are automatically adjusted by periodically sending training 
symbols and allowing the equalizer taps to adjust its parameters in response to these known 
symbols and corresponding received signals. 


Figure 4.9 illustrates how training sequences are organized in a transmission sequence. A 
training sequence with LTg symbols is stored in front of the data package and sent together 
with data repeatedly. Then filter taps are adjusted per iteration by solving the error func- 
tion between these known training symbols and corresponding received signals. A 2x2 SS-SP 
MIMO STE under training mode is shown in Figure 4.10. Where the equalizer 10,1 and w22 
are used to cancel the ISI that is introduced in the direct path channel between Tx;,Rx2 and 
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K{w,, + wA ô (n-L) } 


Filter taps 
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Figure 4.8: Filter taps of a converged equalizer for 2x2 SS-SP MIMO 


Tx2,Rx2 respectively. w12 and wa, cancel the Inter-Stream-Interference between Tx;,Rx2 and 
Tx2,Rx1 respectively. They work together to recover the desired symbols at both receiving 
antennas. 


Total length: L, = Lip + Lp 


Training 
Dat 


Lip Lp 


Figure 4.9: Package structure of trained equalization. 


4.3.1 Zero-Forcing Algorithm (ZF) 
SISO system 


The signal we desire after the equalization operation is the transmitted symbol with an integer 
delay. In other words, after removing both the ISI and Inter-Stream-Interference, an equalized 
signal seg(k) = s(k — ô) where ô is a time delay with an integer number is obtained. Received 
signal can be expressed in a discrete time domain as 


Le 
r(k)= Y h(i)s(k — i) + n(k) (4.8) 


i=— Le 


where s(k) and r(k) are the kth sample of the transmitted- and received signals and 2Le+ 1 is 
the channel length. The ZF (Zero Forcing) algorithm works by forcing the impulse response 


Update 
Weights 


Detection 


Training 
Sequence 


Figure 4.10: The structure for a 2x2 MIMO STE under training mode. 
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of the equalized channel to a delta function. In theory, any channel response can be perfectly 
equalized by a ZF equalizer with infinite number of taps. The upper plot in Figure 4.11 shows 
the channel response through which signals are transmitted. The lower plot is the equalized 
samples which are forced by the equalizer to zeros (2<n<8) except the center tap. From 
Figure 4.4 the equalizer response, due to the channel response H., can be expressed as [20] 


heg = Hew (4.9) 
w = Hjlh, (4.10) 


where he, and w are column vector given by 


w= a (4.11) 
| 


The Channel Response 
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Figure 4.11: The upper plot: Channel response samples. The lower plot: Equalized samples 
using a seven taps ZF algorithm. 


H, is the matrix of channel responses of the form 


he(k) he(k — 1) he(k — 2L) 
NET he(k) he(k — 2L +1) 
H.— | helk+2) — he(k+1) hell — 2L +2) 
| he(k-+2L) h(k-2L—1) e WA | 


Thus 4L samples of the channel response can be used to determine a Zero-Forcing equalizer 
with 2L +1 filter taps. Although the Zero-Forcing algorithm is conceptually simple and 
easy to implement, the impact of the noise enhancement impact is not taken into account. 
Any noise added after the channel gets amplified by a large factor if the channel has a deep 
attenuation and this destroys the overall signal-to-noise ratio. Thus the zero-forcing equalizer 
is not commonly used in most applications. 


4.3.2 Least Square Algorithm (LS) 


MIMO system 


The least square algorithm is widely used to solve over-determined systems, i.e. systems with 
more equations than unknowns. A STE with usage of LS algorithms is for considered and 
equalized signal can be expressed as 


mk) = 9 rj (E) wi (E) (4.12) 


where M denotes the number of receiving antennas, r; is the vector of received symbols at jth 
receiving antenna and Wj, is the filter taps for cancelling the ISI or Inter-stream-Interference 
of the channel hjm. The error function for the mth equalizer output is defined as 


em(k) = 8m(k) — Ym(k) (4.13) 


The index m stands for the mth receiving antenna, (4.12) is repeated for each of Lp outputs 
of mth equalizer corresponding to Lp input signals. Now, collecting Lg samples from the 


mth receiving antenna, which is expressed as rm = [rm(1) rm(2) ... Tm(Lp)]. Then, 
Lg samples of input signals from the nth transmitting antenna is collected and expressed as 
Sn =[sn(1) sn(2) ... sn(Lg)"]. Furthermore, the linear equalizer matrix can be written 
as 
Wii  W241 WM 
Wi2  W22 WM,2 
W = i : (4.14) 
WiM W2M --- WM,M 


where M denotes the number of receiving antennas. In matrix (4.14) the i, jth element of the 

linear equalizer matrix is a linear transversal filter with maximum length of 2L+1 as described 

in Figure 4.4. Thus the minimization problem over a block of Lg samples is expressed as 
Wm = arg min ||Sm — Rw;; ||? (4.15) 


where wm denotes the mth column vector of filter matrix W and || - || is the vector norm 
operation.s,, is the mth column of S matrix as defined in (4.22) and R is introduced below. 
This minimization problem has a unique solution, given by solving (4.16) 


Wm = (RE. R) -- RE. sm (4.16) 


Sm is a vector of source signals of mth transmitting antenna, and R is formed as 


R = [R.,..., Rm] (4.17) 


where, 


Tm(Lt +1) Ta.) << Tm(2) Tm(1) T 0 


um Er Tm(L;—1) ... rat) 0 TE 0 | 
COE T TM MT UT EUN os sat | 
) 


where Tm(k) is the kth sample received at Rx antenna m. 2L4-- 1 is the equalizer filter length 
and L4 is the index of central tap of the equalizer. M is the number of receiving antennas. 
A derivation of the least square algorithm based on the MIMO equalizer can be found in the 
Appendix A. 


4.3.3 Least Mean Square Algorithm (LMS) 
MIMO system 


The LMS (Least Mean Square) algorithm was invented in 1959 by Stanford university profes- 
sor Bernard Widrow and his first Ph.D student, Ted Hoff. It uses stochastic gradient descent 
to start at some arbitrary initial point and iteratively moves towards the optimal point in 
small steps. Due to its simplicity, the LMS algorithm becomes very popular in practical 
applications. Furthermore, it can be applied in both training and DD (Decision-Directed) 
modes. In training mode, as shown in Figure 4.9, by-the-receiver known symbols are added 
and transmitted in front of the data in order to update the equalizer using LMS. In DD 
mode, it assumes that the hard decisions of the received signals are correct and those es- 
timated symbols can be treated as training symbols to update the filter. Using the SGD 
(stochastic gradient descent) the filter updating process is expressed as 


W(k +1) = W(k) — uV.,, Js (4.19) 


where V,J'™® is the gradient of the cost function that is to be minimized, and y is the step- 
size parameter that controls the convergence speed of the LMS algorithm. The remaining 
part to compute is the MSE cost function of LMS algorithm J'™’. The basic idea behind the 
LMS algorithm is to minimize the MSE (Mean Square Error) between the estimated signal 
and desired signal (for simplicity, the system delay ô is ignored ) 


em[(k) = Ym(k) — sn(k) (4.20) 


where n,m stands for the nth transmitting and mth receiving antenna respectively. Here 
n — m is set then the MSE cost function can be written as 


ims E E {lem(k)|?} 
= E (|y (X) — 8n(k)|?} 
= E {rawi = 8n(k)|? } 


where rz is the kth row vector of receiver matrix R which is the same as (4.18) and Wm is mth 
column vector of equalizer matrix defined in (4.14). It is worthwhile to note that the LMS 
algorithm suffer from noise enhancement problems since noise impact is included within the 
criterion so that the minimization of MSE includes the minimization of ISI terms plus noise 
power [43]. The following equations can be derived based on the above J!"* cost function 
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Juas — E {lem(k)|?} 


= E {em(k)em(k)*} 
= E{(rkWm — $n(k)) (teWm — ir D} 
= E{(teWm — sn(k)) (ti Win — sh (k))} 


= E [Ty Wm TE wh Hr wn st (k " — - su (Erw! + su (k)s (k)) 
= wg E Iran wH — wmE (ris? (k)} 
-wm E (ss (k)r jj + E {8n(k)s7(k)} 


E {rar} } is the received signal autocorrelation matrix and E {s,,(k)s*(k)} is desired signal 
variance. Finally, substitution of later term in (4.19) for MSE cost function is done by fol- 
lowing equations 


Val" =E (ril) W-E {s(k)rj} 


W(k +1) = W(k)— rer e W () — s(k)) (4.21) 


where the matrix R and W are defined in (4.18) and (4.14) respectively. s(k) is the kth row 
vector of the matrix S that is defined as 


(4.22) 


Eo "NM KT 


where, N denotes the number of transmitting antennas. 


Stability 

The above equation (4.21) is achieved using LMS (Least Mean Square) algorithm, the step 
size value u has to be chosen so that it satisfies some condition to assure stability of the 
algorithm [21]. This condition is 


O<p<=w (4.23) 


A maz 


where Amaz stands for the largest eigenvalue of the autocorrelation matrix E (nad. p de- 
termines the convergence speed, a large step-size value leads to a fast convergence but high 
residual MSE. On the other hand, a small step-size makes the algorithm converge to relatively 
low MSE but slowly. Therefore, one has to be careful when choosing the step-size, as a wrong 
setting might ruin convergence due to a possible divergence problem. 


4.3.4 Recursive Least Square Algorithm (RLS) 
SISO system 


The Recursive Least Square algorithm, which is similar to the Least Mean Square algorithm, 
aims for minimizing the sum of the squares of the error (the error between the equalized and 
desired signal). By adding a weighting factor, RLS has much faster convergence speed and 
lower MSE compared to LMS [17]. The error between the equalized and desired signal is 
written as 


e(k +1) 2 s(k +1) —w(k)r(k +1) (4.24) 


where the equalizer filter vector w and the received signal vector r have L (filter length) 
elements. s is the transmitted symbol and w(k) denotes the filter taps for the kth iteration. 
Then, the cost function for RLS can be written as 


k 
J'5(k) = V ^ s(t) — wP(k)r(i)]? (4.25) 
i=1 


where 0 < A < 1 is the weighting factor and it is fixed for all k. The cost function is minimized 
by taking the derivatives for all entries of vector w and setting the equation equal to zero [17]. 


k 


= SOs) — w™(k)e(i))r(i) =0 (4.26) 


i=1 


OF™8(k) 
Ow(H) 


In short, w(k) can be achieved and written from above equation as 


w(k) = Ros, (E) ros, (E) (4.27) 


where, 


k 
Reeth) me NO AS EPOD 
i=1 


k 
rcorr(k) == SMA s(i)r(i) 
i=1 


The basic goal of the recursive algorithm is to derive the equalizer filter w(k + 1) based on 
w(k) which leads to the following recursive relations [17] 


Reorr(k + 1) = ARcorr(k) + r(k + 1)rE (k + 1) (4.28) 


and, 


Yeorr(k + 1) = Arcorr(k) + s(k + 1)r(k +1) (4.29) 
From (4.27) it can be obtained that 


w(k +1) 2 Rl, (k + 1)Qrcorr(k) + s(k + 1)r(k + 1)) 


where, 


Arcor (k) = (Reorr(k + 1) — r(k + D)rÉ(k + 1))w(k) 


and, 


w(k +1) = w(k) + Ra. (k + 1)r(k + 1) [(s(k +1) 7 rE(k + 1)w(4%))] 


The matrix R¿,! .(k +1) in the above equation can be updated recursively with the help of 


corr 
the matrix inversion lemma [17]. Given matrices A, U, B and V 


(A+ UBV)7! = A™! — A1U(B71 + VAU) IVA! (4.30) 


By using the matrix inversion lemma, no matrix inverse process is needed in (4.30) which 
reduces the computation and complexity. A, U, B and V are chosen by the following 
equations 


A= ARcorr(k) 
U=r(k+1) 

B-I 

V=rH(k+1) 

then the matrix inverse lemma yields [17] 
1 d H -1 
ET = [Ro Reorr (E) r(E--1)r" (641) Reorr(k) 
Rover (4 + 1) = A [Fez (2) A rH (k+D) Roo, (k)r(k+1) (4.31) 


In order to update the filter taps by following (4.31), an invertible matrix R;,!,.(0) is needed. 


corr 
MIMO system 
In the 2x2 MIMO system, the above equations are modified as 


Rz. Voto (k) 


1 
—1 E —1 
Ross (A T 1) mni A Rs. (Kk) A+rk+1Rcorr(k)ri 4 (4.32) 
W (k +1) = W(k) + Ra (k + Dri [s( +1) — re+1W(k)] (4.33) 


where rẹ stands for the kth row vector of matrix R that has the same format as the matrix 
in (4.18) and W is the equalizer filter taps matrix that has the same format as the matrix in 
(4.14). 


4.4 Blind Equalization 


Sending training symbols at the beginning to initialize the receiver can be impractical and 
not effective for some applications [19]. The correct sampling point has to be known in order 
to extract the training symbols before the equalizer converges. However, trained algorithms 
need information from training sequences to converge. In addition to the above issues, the 
training sequence will also consume bandwidth. So, in some cases, it is desirable to equalize 
the channel without the aid of training sequence and it is called blind equalization ([51], [13], 
[28], [38], [46]). Figure 4.12 illustrates the structure of a 2x2 MIMO equalizer using blind 
algorithms. 
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Figure 4.12: MIMO Equalizer using blind algorithm. 


4.4.1 Constant Modulus Algorithms (CMA) 
SISO system 


Many modulation format have the constant modulus property such as PSK, FSK. The received 
signal r(k) may then be described by 


Lc 
r(k) = M h(i)s(k — i) + n(k) (4.34) 
i=1 


where s(k) is the transmitted signal and r(k) is the received distorted signal plus noise and Le 
is the channel length. The source is unknown but it has a constant amplitude value |s| = C 
for all k. So, the objective is to construct a receiver weight matrix w such that equalized 
signal y(k) can be denoted by 


y(k) = w™(k)r(k) (4.35) 


where w(k) is defined as w(k) = [w(k — 5), w(k—£+1),...w(k), w(k+1),....w(k+ UN and L 
is the filter length. r(k) is defined as r(k) = [r(k— 3), r(k— £41), nmr(k), r(k4-1), cr)" 
Thus the possible solution is to compute a w such that |y(k)| =C for all k where y(k) is 
the output of equalizer. Again, a stochastic gradient descent (SGD) method is considered 


w(k+1) = w (k) —uVyJ™ (w) (4.36) 


y= ig { Com Ros) } (4.37) 


where Rema is known as Godard radius [43] and it is a real constant and the value is cho- 
sen based on source alphabet and p factor. It was noted [43] that by setting Rema = 
E {|s(k)|??} /E {|s(&)|?} and p = 2 the existence of a local minima of VwJ°™ that per- 
fectly equalizes the distorted signals with an arbitrary phase rotation is assured in a SISO 
system. By computing the gradient of the cost function we obtain (4.39) for updating the 
weights [38] [41] 


Vw Som = y(k)(ly(k)|? — Rema)" (k) (4.38) 


w (k + 1) = w (k) — uy) (ly) — Rea)r" (k) (4.39) 


MIMO system 


In the MIMO case filter taps are updated independently based on different equalized signals 
from each receiving antennas 


Ym(k) = rxwm(k) (4.40) 


Wm(k + 1) = Wm(k) — paw) ym (ke) |? = Roma) Ek (4.41) 


where w,,(k + 1) is the mth column vector of the equalizer matrix that is defined in (4.17), 
and Ym(k) is the mth equalized signal for mth receiving antenna. rp is the kth row vector of 
receiving matrix R that is defined in (4.18). It is noteworthy that the CMA not only removes 
ISI but is phase blind. So, the equalized signals after CMA equalizer might have a constant 
phase rotation that needs to be rotated back. Additionally, as shown in [11], the CM criterion 
of (4.37) also works well for non-constant modulus cases such as 16QAM or 64QAM. But the 
minimization of (4.37) does not ensure the recovery of all source signals in a MIMO system 
because it may converge to recover the same source signal at many outputs. In order to solve 
this problem, a cross-correlation term is introduced due to its computational simplicity [32]. 
Then (4.38) is rewritten as 


m-—1 
Vw, J^" (wm) = (1 gm (E) (lys (E). — R)rg +a 9 milk) yl) re, (m = 1,2, ...M) 
i=1 


(4.42) 
equation (4.39) will be modified as 
m-—1 
Win(k + 1) = Win(k) — uy (E) — 0) (lum)? — R) +a 7 Fmi(kyi(kIre — (443) 
i=1 
Pmilk +1) = APmi(k) + (1 — A) Ym (K)y; (k) (4.44) 


where a is the mixing parameter which is chosen between 0 and 1. £,;(k) = E lym(k)y;(k)] 
is the cross-correlation between the mth and the ith equalizer outputs and prevents the 
extraction of the same signal at many outputs. A € [0,1] is a parameter that controls the 
length of the effective data window in the estimation. The first term in (4.42) ensures the 
recovery of only one signal at each equalizer output and the cross-correlation term ensures 
that each equalizer output is different from other ones. œ has to be selected carefully for 
its optimum value in different channel conditions. For example, higher o value is chosen 
if equalizers are more likely to converge the same channel. Here if we assume a 2x2 Los 
MIMO system with an optimal separation (Inter-Stream-Interference from Tx has a 90° 
phase shift compared with desired transmitted signals at Tx). Initial filter taps are set up 
by the following matrix 


0 0 
Wi1 e? E 021 
0 0 
DE =) a 
| ae | (4.45) 
ii ji 
EN 
w12 e$ e | W22 
0 0 
blo 3 l 


where wy = [wi, Wi]? and wz = [w»21 wa9]'. Each column vector wij has a length of 
L taps. el, e13 e Iz ande? are the central taps for w11, w12, wo1andw»» respectively. There 
is another algorithm that uses the Gram-Schmidt orthogonalization procedure to force this 
separation, more detailed discussion can be found in paper [31] for reader's interest. Figure 
4.13 shows the initial filters setting in the simulation system (also given in 4.45), where it is 
assumed that the 2x2 SS-SP MIMO system starts with an optimal antenna separation. 


60 


K{w,, +w,, 8 (n-L)] 


Real magnitude 
o 
a 
T 


o 


Filter coefficents (43taps) 
S {w *WA 6 (n-L) } 


o 


Filter coefficents (43taps) 
R {w,, + Wp Ò (n-L) } 


Imaginary magnitude 
5 
a 
T 


Real magnitude 
o 
a 
T 


o 


Filter coefficents (43taps) 
3 (w,, + W,, 6 (n-L) } 


o 


Imaginary magnitude 
5 
a 
T 


Filter coefficents (43taps) 


Figure 4.13: Filter initialization of STE for the 2 x 2 SS-SP MIMO 


4.4.2 Modified Constant Modulus Algorithm (MCMA) 


SISO system 

The CMA (Constant Modulus Algorithm) algorithm expects constellations with a constant 
amplitude, which degrades the equalization performance when a M-QAM (M>4) modulation 
is used since it has multi-modulus property. In order to improve the performance of the CMA 
for M-QAM (M>4), a multi-modulus algorithm, called MCMA (Modified Constant Modulus 
Algorithm), has been proposed in [29]. In this algorithm the cost function for CMA (4.38) is 
modified to consider the real and imaginary parts separately. The modified cost function is 


written as 
Je? k) = Jrk) + Jr(k) (4.46) 


where Jr(k) and Jr(k) are the cost functions for the real- and imaginary parts of the equalizer 
output y(k) = yr(k) + j - yr(k) respectively and they are defined as 


In(k) = EQ GG) — Re)’ } (4.47) 
40) = E{(\S{y(e)} ? - my) (4.48) 


Assuming the input data is complex numbers with two dimensions, Rp and R are the real 
constants determined for the real and imaginary parts of the source signals respectively 
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ICONO 

Er = EQBHSQOE) (4:49) 
ECONO: 

Br = ASEO Eu 


where 9e (s(k)) and S (s(k)) denote the real and imaginary part of source signals s(k) re- 
spectively. The tap weights vector is updated using SGD to optimize the cost function 


w(k+1) = w(k)— u Vw J^"? (w) (4.51) 
= w(k) — ue(k)r (k) (4.52) 


where the error signal e(k) = er(k) + j * e7(k) is given by 


en(k) = R{y(k)} (R {y(k)} P — Rr) 
e(k) = SQ) S {y(k)} |? — Ri) 


MIMO system 


For MIMO systems, equalizer taps are updated independently for different receiving anten- 
nas which is similar to the Simplified Constant Modulus Algorithm (SCMA) method. This 
processing is expressed as 


Ym(k) = PkWm(k) (4.53) 


emR(k) = R{Ym(k)} (IR {ym(k)} |? — Rr) (4.54) 


( 
( 


milk) = S{Ym(k)} (IS {Ym(k)} |? — Ri) (4.55) 
em(k) = emprl(k) +j *enalh) (4.56) 
Win(k +1) 2 Wm(k) — em(k)rÉ (4.57) 


here rz is the kth row vector of receiving matrix R that is defined in (4.18) and ym(k) is 
the kth sample of equalized signals on the mth receiving antenna. This algorithm modifies 
the CMA cost function so that the resulting cost functions becomes dependent of the phase 
rotation. Hence, the MCMA can correct for some carrier frequency offset and remove the 
ISI simultaneously. As a result MCMA results in a performance improvement in convergence 
speed and residual MSE floor [29]. 
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4.4.3 Simplified Constant Modulus Algorithm (SCMA) 
MIMO system 


A Simplified Constant Modulus Algorithm (SCMA) has recently been proposed in [32]. In- 
stead of projecting signal points onto a M-QAM constellation based on both real and imagi- 
nary part, it only considers one dimension (either real or imaginary part due to their symme- 
try) which is more flexible and less complicated compared to CMA or MCMA from practical 
implementation point of view. This algorithm suggests to minimize the following cost function 
for the mth Rx output and based on the real-part of signal in this thesis 


Ja =E (C {Ym(k) } |? E Rima)? (4.58) 
where R{-} stands for the real-part operation. R is given as 
p ERGO 

E {IR (50) |?} 


In order to force the separation of equalized signals at different outputs, the same cross- 
correlation term is used here as in (4.42) and the gradient of (4.58) is rewritten as 


(4.59) 


Yml k) = rkWm (Kk) 


m-—1 
VRE = (L—a)R {ym(k)} UR (o (K)) |? — R) +00) 5 Pmalk)yi(k) bre 
i=1 
m-—1 
Winlk +1) = Wilk) — ut fy (k)} A — UR Qu (k)} |? — R) +a V 5 Ami (k)yi(h) Ire 
i=1 


where a is the mixing parameter which is chosen between 0 and 1. And fmi(k) = E [ys (k)y; (k)] 
is the cross-correlation between the mth and the ¿th equalizer outputs and prevents the ex- 
traction of the same signal at many outputs. Ym(k) is the kth equalized signal of mth receiving 
antenna. rz is the kth row vector of receiving matrix R that is defined in (4.18). The cri- 
terion of SCMA leads to the converge of the signal constellation that, in absence of noise, 
corresponds to recovery of desired signal with a constant phase rotation that differs between 
modulation formats. For QAM this rotation is 5. It is shown in [32] that the proposed 
algorithm presents a lower computational complexity compared to the MCMA without any 
loss in performance but it can be less stable than the other two algorithms in some specific 


channel conditions. 


4.5 Decision-Directed Mode (DD) 


In addition to blind and training mode, there is another mode called Decision-Directed mode 
that is commonly used. The DD (Decision-Directed) mode that has been implemented for a 
2x2 MIMO system is illustrated in Figure 4.14. The structure only shows the output of the 
equalizer for the first receiving antenna. The main idea of DD mode is that it assumes the 
current hard decision of equalized symbol from previous iteration is correct and the channel 


63 


is not changed during this iteration period so that filter taps from previous iteration can be 
used to equalize the received signals for the current iteration and treating equalized symbols 
as training symbols to update the filter taps for the current iteration. Then the training 
mode algorithms such as LMS and RLS processes equalization without receiving any training 
sequence from the transmitter. This is why blind equalization is normally followed by a DD 
mode in order to significantly improve the MSE of the estimation. Thus in DD mode perfect 
initial convergence is not needed for the equalizer to converge to the optimal setting as long 
as estimated signals are correct most of the time. 
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1 " 
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Y Y Y Y y (k) 
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Figure 4.14: 2x2 MIMO-STE under Decision-Directed Mode. 


MIMO system 


For RLS algorithm the error function of (4.24), in DD mode, can be rewritten as 


emlk + 1) = 8m(k + 1) — rr+1Wm(k) 


where rz41 is the k + 1th row vector of R that defined in (4.18). Wm(k) is the mth column 
vector of W matrix that is defined in (4.14). Sm(k + 1) is the hard decision of equalized 
signal y (k +1) at sample time k+ 1 on the mth receiving antenna. Then, filter taps can be 
updated by following equation 


Wk +1) = Wk) + Ras, (e + Urey BE +1) — ri W(4)] 
For LMS algorithm, the error function is modified as 


emlk) = 8m(k) — rkWm (k) 


where ŝm(k) is the hard decision of equalized signal ym(k) at sample time k for the mth 
receiving antenna. Equation (4.21) is rewritten as 
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W(k +1) = W(k) — star e W(K) E 


4.6 Complexity comparison of different algorithms 


Since there is always a trade off between the performance of the algorithms and their compu- 
tation complexity. It is important and necessary to make a complexity comparison between 
different algorithms [31] and try to find a algorithm that might achieve the best balance. 


Algorithm | Multiplications Additions 

CMA 2(4M-L+3)-M (9M-L-2). M 

MCMA 2(4M-L+3)-M (9M-L+5)-M 

SCMA (1AM-L+3)- M (6M-L+1) - M 

LMS (SM -L+2) - M (9M-L-2). M 

RLS 4M?.L(1+6L+M-L) | 22M? .1?+2M+1+10M?.L+4M-L 


Table 4.1: Comparison of complexity of different algorithms per weight update. 


In Table (4.1), L denotes the order of equalizers (number of filter taps) and M is the number 
of receiving antennas. From this table it can be seen that SCMA has the lowest computational 
complexity. It is noteworthy that (as shown in the simulation section), SCMA has almost 
the same performance as the MCMA algorithm but lower computation complexity. This 
might be a big benefit for the real implementation since SCMA only use the information from 
one-dimension of signals to update weights instead of two. Compared with LMS, the RLS 
algorithm has a faster convergence speed but much higher complexity which is shown in the 
above table. 
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Chapter 5 


Simulations 


The simulation software was written in Matlab using a modular architecture to enable easy 
selection and combination of different blocks. A graphical user interface (GUI) was con- 
structed alongside the code to ease block selection and simulation parameters selection. The 
simulator handles both 2x2 SS-SP LoS MIMO and 4x4 SS-DP LoS MIMO. All the simula- 
tions are done at baseband using Matlabs 64 bit floating point numbers making quantization 
errors negligible. A block diagram of the entire 2x2 SS-SP LoS MIMO system can be seen 
in Figure 5.1. It gives a good overview of the system design and makes it easier to visualize 
how data is processed in the system. Both branches generate independent data which remain 
separate until being combined by the MIMO channel. At the Rx the data streams are pro- 
cessed separately with the exception of the STE which needs both data streams in order to 
cancel the inter-stream interferences. Each Rx branch then outputs its estimate of the trans- 
mitted symbol- and bit sequence and thus the branches have their own BER and MSE values. 


Following the sequence of events necessary to generate and process a block of data bits gives 
a better insight into the design. Data bits are generated using a uniform random distribu- 
tion which are modulated in the modulator into a Gray coded M-QAM constellation. The 
resulting symbol sequence is then filtered using a oversampled root-raised-cosine filter. The 
filter emulates an analog signal by an oversampling rate of 16, resulting in a symbol time 
of Tos = (Q- Fs)? = 1.25 ns, and a roll-off factor of 8 = 0.35. The signal then passes 
through the LoS MIMO channel which will be discussed in Section 5.3. At the receiver ad- 
ditive white Gaussian noise with zero mean is added to the signal. Independent realizations 
of the AWGN is added to each signal stream. This is then followed by a matched filter that 
is not oversampled as in the Tx. Phase noise with a certain spectral density Spp(w) is then 
generated once (new realization for each symbol) and multiplied to both signal streams which 
is consistent with the assumption of a synchronous system (one phase noise realization for the 
whole system). The analog to digital converter then samples the signal at baud rate without 
any timing errors. The signal then enters the space time equalizer (STE) where it is equalized 
and the inter-stream interferences are cancelled. The STE tunes its equalizers by using the 
error signals generated by the detectors. The detectors have phase trackers that rotate the 
signal so to alleviate some of the phase noise. The sliced signal is then finally feed to the 
demodulator which outputs the estimated data sequence. 


The following approximations where made for the design of the system (most of which can 
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easily be dropped in future versions). 


e The system is considered to be completely synchronous so that there is no need for a 
symbol timing estimation at the ADC and no need for any phase adjustments between 
the Rx signal streams in the STE. There is also no carrier frequency offsets between Tx 
and Rx. 


e There are no quantization errors introduced to the system. The simulation is done at 
a 64 bit floating point number representation making it highly accurate. 


e Between every Tx and Rx pair there are only two paths considered, one direct LoS path 
and one ground reflection. The direct paths all have the Rx power of one and all the 
ground reflected paths have the same amplitude b. All the ground reflections have the 
same delay as well 7 = 6.3 ns. This is the Two Ray model. 


e The Tx and Rx antennas can only move back and forth in a plane that is parallel to 
the masts. 


These assumptions ease the design of the simulation and allow for more easily interpretable 
results. 


STE Detector 


Figure 5.1: A block diagram representing the tentative design of the simulator for a 2x2 
SS-SP LoS MIMO. 


A 4x4 SS-DP LoS MIMO system will also be investigated. The design is almost the same as 
shown in Figure 5.1 with the exception that there are four branches now at the Tx and Rx 
and the design of the STE now includes 16 equalizers (including 8 XPICs). This is shown for 
one polarization at the Rx in Figure 5.2. The dual polarized 4x4 channel matrix is modeled 
using (3.25) which assumes that both the Tx and Rx have two dual polarized antennas (both 
polarized antennas collocated). Each ground reflection has though an independent phase off- 
set creating 16 independent frequency selective channels. 


For convenience some basic simulation constants are summarized in Table 5.2. All of them 
are kept constant during the simulations unless otherwise stated. 
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Figure 5.2: A block diagram representing the tentative design of the vertically polarized Rx 
for a 4x4 SS-DP LoS MIMO. The block diagram neglects all noise sources. 


5.1 Goals 


The goal of the simulation is to evaluate the performance of a LoS MIMO system, both 2x2 SS- 
SP and 4x4 SS-DP, subjected to different conditions and using different adaptive algorithms. 
The results should then work as a guideline for an implementation of a real system and help 
to select algorithms and give an understanding of the impact of certain suboptimalities in 
the installation of the system as well as response to different channel conditions. To be more 
specific the desired results are given in Table 5.1. 


Section | Description 

6.2 Tolerance to low XPD 

6.2 Tolerance to suboptimal antenna placements 
6.6 Tolerance to phase noise 

6.1 Tolerance to frequency selective channels 

6.5 Tolerance to mast swing 

6.4, 6.3 | Evaluation of different algorithms for STE 


Table 5.1: Desired results obtained by simulation and the section where they are published 
in the report. 


The evaluation of different adaptive algorithms (STE's) is perhaps the biggest task. For that 
part the main focus will be on blind equalizers as well as training based methods to initialize 
the system followed by Decision Directed (DD) adaptive algorithms to track the channel. For 
these methods there are mainly three attributes that are investigated 


e Convergence speed 
e MSE- and/or BER-floors 


e Tolerance to mast swing and phase noise 
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5.2 Quantifying Performance 


Comparing performance requires some measurable quantities that can be used to evaluate 
different setups and different channel conditions. There will be mainly three quantities used 
for comparison. The first one is the bit-error-rate (BER) defined as the ratio between the 
number of bits received in error to the total number of bits received or more precisely 


Ly 


where bn is the nth transmitted bit (0 or 1) and bn is the nth estimated bit at the receiver. 
Ly is the total length of the bit sequence. This is a very practical quantity to measure for 
a real system since that will define for a large part the actual performance of a real system. 
Since the simulator does not use any error correction coding the BER can easily be lowered 
by using a good coding scheme. The target BER is though often very low requiring very long 
simulation times to get obtain reliable estimates. Another quantity that does not require such 
long simulation times is the mean-square-error (MSE). It averages the square error between 
the soft estimate of a symbol and the actual symbol transmitted 


SLE gn fal 


B 


MSE — 


Where sn is the nth transmitted symbol and f; is the soft estimate of the nth received sym- 
bol. Lg is then the number of symbols transmitted. 


The final quantity is the penalty described in Section 3.6. It quantifies how much more the 
Tx power needs to be increased between different scenarios to achieve the same BER and is 
usually specified in dB. The target BER for the penalty will be 107%. A baseline SNR needs 
to be established for the penalty, that is usually the lowest achievable SNR for a certain BER. 
For a SISO system with a flat channel and using 16QAM modulation the SNR is about 10.5 
dB for a BER of 1073. A 2x2 MIMO system will have a 3 dB array gain! so the baseline 
SNR will be considered to be 7.5 dB. Finding the penalty using simulation software can be 
a challenge since the BER is an output and is a function of great many variables. The SNR 
needed for a certain BER is desired to be found for a certain setup. For these simulations a 
search algorithm inputs the SNR into the simulator and runs it successively until the desired 
BER has been found. That SNR minus the baseline SNR is then the penalty. 


5.3 The Channel 


The frequency selective channels used by the simulator is the Two Ray channel described 
in Section 3.4. The Tx and Rx antenna discs are mounted on two masts of height Lr and 
Lg respectively. The Rx mast is situated on a hill of height Lz and the horizontal distance 
between the base of the masts is Ly. The masts can then have a time varying tilt of T for 
the Tx and fig for the Rx to simulate mast movement in heavy wind. The setup can be seen 
in Figure 5.3. Unless otherwise specified the parameters for the channel setup are the same 
as given in Table 5.2. 


Two Rx antennas doubles the received power. 
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Figure 5.3: The setup of the Tx and Rx antennas. Pr, Br, dr, dr, Lr, Lg, Lr and Ly can 
all be varied. 


The amplitude of the reflected waves is b and for the purpose of the simulation it is bounded 
between zero and one where one is the amplitude of the direct wave. Each ground reflection 
must traverse a longer path so there is a delay difference between the ground reflection and 
the direct path. This delay is set to be fixed in accordance to Rummler’s model at 6.3 ns 
which is a common standard but the phase offset of the arriving direct path and ground 
reflection is random. For a 2x2 MIMO system there are 4 channels where each one has an 
independent phase offset between ground reflection and direct path but otherwise share the 
same parameters. For example the channel between Tx2 and Rx, will have the same delay 
and same ground reflection amplitude as any other channel but it has an independent phase 
offset $521. This phase offset is enough to determine how the channel behaves. It can be 
relatively flat or it can act as a lowpass-, highpass- or bandpass filter. If the channels have a 
deep notch then noise will be amplified but attenuating channels do not necessarily represent 
the worst case scenario. Since the channels produce collectively some frequency dependent 
MIMO channel with some frequency dependent singular values as is discussed in 3.3 which 
influence the performance it is not always easy to identify the worst case channel. 


Since there are infinitely many channel realizations it is necessary to narrow down the se- 
lection to a few cases. Selecting the most interesting channels can be challenging since it is 
not even clear what represents a worst case. The selection process thus becomes somewhat 
ad-hoc. The channel responses displayed for the MIMO channels are obtained by measuring 
the impulse responses seen between each Tx- and Rx antenna pair. This is done by trans- 
mitting an impulse on one transmitter at a time instant and measuring the output from each 
receiver. The impulse responses include the Tx- and Rx filters as well as the ADC. Selecting 
the channels for the 4x4 MIMO is clearly even harder since there are 16 independent chan- 
nel responses making up the whole MIMO channel. To ease the selection the 4x4 MIMO 
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channels are constructed from the 2x2 channels making the following simplifications that 
both the polarizations and corresponding cross-polarized paths see the same channels as the 
single polarized case. That is the channels between Tx, and Rxm are the same for the SP 
and DP cases independent of polarization (hnymy = hagag = ahnymy = ahnymy where 
a is some constant related to the XPD). This is made more apparent in the following sec- 
tions where the two channels chosen for use in the simulations are described (frequency plots 
are obtained using DTFT). They are denoted A and B and both have their defining attributes. 


These channels are also realized as SISO channels by removing the paths h,m where n £ m. 
The BER for the resulting parallel channels are then averaged for penalty calculations. 


5.3.1 Channel A: All Notches 


2x2 SS-SP LoS MIMO: This channel has a notch between all Tx- and Rx pairs. This will 
attenuate part of the spectrum of the signal by the full notch depth with a subsequent SNR 
reduction at the receiver. The placement of the notches can be seen in Figure 5.4(a) where 
the amplitude response is plotted from —F;/2 to Fs/2 where Fs is the sampling frequency. 
The depth of these notches depend on the amplitude of the ground reflection selected here to 
be b — 0.9. The channels act as FIR filters and have thus a guaranteed linear phase response 
which can be seen in Figure 5.4(b). Looking at the pole zero plot for the z-transform of the 
impulse responses in Figure 5.4(c) it is noteworthy that all of the channels have a zero rather 
close to the unit circle making them harder to invert at the Rx. Finally the condition number 
along with the singular values is plotted in Figure 5.4(d) as a function of frequency. 


4x4 SS-DP LoS MIMO: This channel is now significantly more complex even though 
it is constructed from the channel in Figure 5.4. The magnitude of the impulse response is 
plotted in Figure 5.5 while the condition number and singular values are plotted in Figure 
5.6. 
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Figure 5.4: Channel A: Frequency response of the four SISO channels making up the 2x2 
SS-SP LoS MIMO channel along with pole zero plots and a plot of the condition number and 
singular values as a function of frequency. For these plots b — 0.9. 
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Figure 5.5: Channel A: The magnitude of the impulse 
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Figure 5.6: Channel A: The condition number and 
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responses for the 4x4 SS-DP LoS 


singular values for the 4x4 SS-DP LoS 


5.3.2 Channel B: Ill conditioned 


2x2 SS-SP LoS MIMO: This channel has only one deep notch but is very ill conditioned 
since the condition number becomes very high for some part of the spectrum as is shown in 
Figure 5.7(d). The magnitude response is seen to contain only one deep notch while the other 
channels are relatively flat (Figure 5.7(a)). The phase is again linear as is expected and as 
can be seen from Figure 5.7(b). Finally the pole zero plot shows that the channel containing 
the notch has a zero closest to the unit circle (Figure 5.7(c)). 
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Figure 5.7: Channel B: Frequency response of the four SISO channels making up the 2x2 
SS-SP LoS MIMO channel along with pole zero plots and a plot of the condition number and 
singular values as a function of frequency. For these plots b = 0.9. 


4x4 SS-DP LoS MIMO: The magnitude of the impulse responses is plotted in Figure 
5.8 while the condition number and singular values are plotted in Figure 5.9. Some of the 


notches in Figure 5.8 seem to be deeper than others, this is because of the XPD chosen to be 
10 dB for the plots. 
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Figure 5.8: Channel B: The magnitude of the impulse responses for the 4x4 SS-DP LoS 
MIMO channel. Here XPD is 10 dB. 
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Figure 5.9: Channel B: The condition number and singular values for the 4x4 SS-DP LoS 
MIMO channel. 
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5.3.8 Mast Swing 


Mast swing will cause the antennas to move and this movement will change the channel. The 
length of the direct path will change which in turn changes the phase of the direct path. 
Since the ground reflection phase is modeled as static this will change the offset between the 
phase of the ground reflection path and the direct path causing the notches seen in Figure 
3.11 to move across the channel. This forces the STE to adapt continuously to new channel 
conditions and this movement sets a lower limit for the convergence speed of the adaptive 
filters. 


How the channel behaves with mast swing in reality is poorly understood since the phase 
offset between the LoS path and the ground reflection path evolves with an unknown distri- 
bution as the antennas moves. If the model was entirely calculated from the geometry then 
the conclusion could easily be drawn that the phase offset hardly changes at all. A millimeter 
wave reflecting off some unknown surface will though experience a phase shift that can differ 
even for only slight changes. This is enough to motivate that it is not prudent to assume 
a static channel when the antennas are moving. Fixing the phase of the ground reflection 
while varying the LoS paths according to geometric calculations seems to be an intuitive 
compromise representing a "bad enough” scenario. This is further motivated by the fact that 
the ground reflected path is modeled as having a static relative delay. Quantifying how much 
a standard mast might move in heavy wind can also be quite difficult. The values chosen in 
table 5.2 are only motivated compromises explained below. 


There have been measurements of the width (in degrees) of half of the 3-dB beam widths for 
parabolic antennas commonly used for microwave links. There is quite a selection of antennas 
with different gains intended for different frequency bands. A high half 3-dB beam width is 
here considered to be 2.3°. The value 2.3° is thus chosen to represent the highest possible tilt 
of a mast. It is a well known fact that the frequency of oscillation for a pendulum made up 


of a solid uniform rod of length L is 
1 3g 
f= on 2L 


where g is the acceleration of gravity (9.8 [m/s?]). For a 10 m long rod the result is f ~ 0.2 
Hz. This value is used as the maximum frequency of mast swing for the simulations. 


5.4 Space Time Equalizers 


A space time equalizer is implemented at the receiver in order to compensate for inter- 
symbol-interference, inter-stream-interference, cross-polarization-distortion and also to track 
the channel changes. One scope of this master thesis is to analyze and evaluate different 
adaptive algorithms for space time equalizers. Due to the complexity of implementation is- 
sues, only MMSE (Minimum mean square error) algorithms are analyzed and evaluated in 
this thesis. 


According to the properties of the chosen channels, a discrete-time baseband FIR equalizer 
with 43 taps is set up in all simulations and the 21st tap is set as the center tap. 
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5.4.1 Trained Decision Directed Equalization 


The scheme used by the space time equalizers fall into two categories: Blind followed by 
Decision Directed mode and Trained followed by Decision Directed mode. The trained al- 
gorithms are mainly LS (Least Square) but LMS (Least Mean Square) and RLS (Recursive 
Least Square) can also be applied in training mode. The Decision-Directed algorithms in- 
clude LMS (Least Mean Square) and RLS (Recursive Least Square) along with some variants. 


The trained algorithms use 200 and 400 training symbols for 2x2 SS-SP and 4x4 SS-DP MIMO 
systems respectively to obtain the initial convergence. The weights update of equalizers will 
then be switched to DD (Decision Directed) mode. The training sequences are modulated 
from the same pseudorandom binary sequence (PRBS) so that each realization of the filter 
taps generated by training and the same algorithm should be similar (noise will also affect 
the filter). 


5.4.2 Blind Equalization 


Three different blind algorithms will be compared and they are CMA (Constant Modulus Al- 
gorithm), MCMA (Modified Constant Modulus Algorithm), and SCMA (Simplified Constant 
Modulus Algorithm). Since the blind algorithms are phase independent it is assumed that 
the blind algorithms are accompanied by a phase-recovery loop that rotates the constellation 
into its correct place and this rotation acts per 100 symbols. It is also assumed that there is 
no carrier frequency offset between the transmitter and the receiver. After the initial conver- 
gence of the blind algorithms, the STE will switch into Dicision-Directed mode to be able to 
track the channel changes and improve the MSE. 

The initial filter taps are set based on the assumption of optimal antenna separation which 
leads the inter stream interference from another transmitting antenna (2 x 2 MIMO) to have 
a 90° phase shift compared to the desired signal. The initialization of the filters are based on 
this assumption. Thus all the central taps of the equalizers will be set to one or ein 
shown in Figure 5.10. 


as is 


5.5 Noise and Phase Tracking 


White Gaussian noise is added to each Rx signal stream separately with a relative power 
measured from one LoS signal stream. Thus the actual noise power is independent of ground 
reflections and how the LoS signals add at each Rx (if it is destructive or constructive ad- 
dition). This results in a 3 dB array gain for a 2x2 MIMO system compared to a SISO system. 


Phase noise is generated in the time domain and is applied on a symbol-to-symbol basis. The 
system is considered to be synchronous and driven only by one oscillator at the Rx. Thus 
only one realization of the phase noise is generated and applied to each stream. The streams 
will thus have fully correlated phase noises. The Welch power spectrum of the phase noise 
is shown in Figure 5.11. It is a 1/f? noise with -85 dBrad?/Hz? power at 100 kHz from the 


? B dB[rad?]/Hz = B — 3 dBc/Hz 
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Figure 5.10: Filter initialization of STE for the 2 x 2 SS-SP MIMO. 
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Figure 5.11: The Welch power spectrum for the phase noise. 


The phase tracker used to compensate for phase noise is shown in Figure 5.12. It is applied for 


78 


each Rx branch as is shown in Figure 5.1. For handling differential phase noise (independent 
phase noise realization for each Rx branch)? NM phase trackers and phase rotators are 
needed. One in front of every equalizer, inter-stream canceler and XPIC. This is due to the 
relatively low bandwidth of these filters. 


I(s) 


Koi Kio 


Figure 5.12: Schematic view of the phase tracker. 


The gains of the phase tracker were found by brute force search. The actual values chosen 
can be seen in Table 5.2. 


5.6 Simulation Parameters 


The parameters used for the simulations are summarized in Table 5.2. They will be kept 
constant during all simulations unless otherwise stated. 


?In a microwave link with radios separated several meters the phase noise will be differential. 
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Symbol | Value Description 

F; 50 MHz The sampling frequency 

Q 16 Oversampling rate at the Tx filter 

Fe 17 GHz Carrier frequency 

B 0.35 Roll-off-factor for the root-raised Cosine filters 

M 16 Modulation order for M-QAM 

Lr 10 m Tx antenna mast height 

Lg 10m Rx antenna mast height 

Ly 1 km Tx and Rx antenna mast horizontal distance 

dr 2.97 m Tx antenna separation 

dR 2.97 m Rx antenna separation 

Lr 0m Rx landscape height 

Bmax 2.3? Maximum tilt of antenna masts 

FMs 0.2 Hz Frequency of mast swing 

L 43 Number of filter taps for each filter in the STE 

n 21 The placement of the central tap in each filter in the STE 
Kn 0.35 A gain for the phase tracker 

Kro ios A gain for the phase tracker 

Kp 0.015 A gain for the phase tracker 

Homa, 5e-6 The stepsize of the CMA. 

Lmcma | 5e-6 The stepsize of the MCMA 

Hims 1e-4 The stepsize of the LMS in the Decision-Directed Mode 
a 1 The mixing parameter for SCMA 

b 0.9 The ground reflection factor 

SNR 35dB Signal-to-noise ratio 

Lir25 200 symbols | Number of training symbols for 2x2 MIMO system 
Ltras 400 symbols | Number of training symbols for 4x4 MIMO system 

R 13.2 The average constant magnitude of Blind algorithm (16QAM) 
W 0.98 Forgetting factor for RLS algorithm 

BER; 1078 Target BER for penalty calculations 


Table 5.2: Parameters used during the simulations. 
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Chapter 6 


Results 


Here the simulation results will be presented. When a relevant value used for a particular 
simulation is not stated in the text the reader is referred to Table 5.2 for a list of default 
values. 


6.1 Frequency Selective Channels 


Two MIMO channels are selected for the simulations, those are discussed in Section 5.3. It 
is interesting to see how the system responds to these channels and try to get a feeling for 
what requirements they put on the design. The length of the STE equalizer filters will heavily 
affect the performance since a too short filter will lead to high- or even infinite penalty. 
Figure 6.1(a) plots the penalty versus filter length for a 2x2 SS-SP LoS MIMO while Figure 
6.1(b) does this for a 4x4 SS-DP LoS MIMO with three different values of XPD. The channels 
were static during the simulations and a long (600 and 2400 symbols for the 2x2 and 4x4 sys- 
tems respectively) training sequence was used to initialize the STE using the LS algorithm, 
the filter taps were kept static after the training sequence. Looking at the plots it can be seen 
that at a certain number of filter taps a penalty floor is reached. The penalty then slowly 
grows as the number of filter taps are increased since each tap introduces a slight estimation 
error. For channel A the penalty floor is reached at about 15 taps while for channel B it 
takes about 43 taps. The penalty floors for channel A and B are then at about 14- and 19 
dB respectively. For a low XPD value it is even beneficial to have more taps but for the 
sake of complexity the number is chosen to be 43 for the simulations. Removing the ground 
reflection completely (b = 0) a penalty of 0.8 dB for the 2x2 MIMO is still observed. This 
is due to the fact that the taps of the STE introduce some estimation errors which add up. 
Increasing the number of training symbols to 6000 reduces this penalty to 0.2 dB. 


It is interesting to see what kind of penalties are expected for a random channel realization, 
that is random phase offsets between the direct paths and ground reflections while the b-values 
and delays are kept fixed. This is done in Figures 6.2(a) and 6.2(b) for 50 random channel 
realizations for 2x2- and 4x4 MIMO respectively. 

The penalty can be negative since the ground reflection can enforce the direct path which is 
seen in some cases. None of the channel realizations exceed the penalty of channel B. 
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Penalty vs. Filter Length. 2x2 SS-SP LoS MIMO 
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Figure 6.1: Penalty (at BER of 107°) as a function of filter length for the STE. 


Figure 6.3 shows how the penalty changes as a function of the ground reflection amplitude 
for both channel A and B. Channel A behaves as might be expected but channel B reaches 
its highest penalty at b = 0.89 and then the penalty reduces as b increases. An explanation 
for this is that the frequency dependent condition number has a peak at b = 0.89. Increasing 
the strength of the ground reflection further starts to improve the condition number. 

It is interesting to compare these results to a SISO system. By removing the diagonal paths 
of the 2x2 MIMO channels (that is removing h1» and h21) two parallel SISO channels are 
obtained. The BER is then calculated as the mean of the two streams. Penalty is then plotted 
against the number of filter taps in Figure 6.4. For 13 filter taps the penalty is seen to be 9- 
and 8 dB for channel A and B respectively. This is a much lower penalty than seen for the 


MIMO channels and with the added benefit of smaller equalizer filters. 
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Penalty for random channels Penalty for random channels. 4x4 
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Figure 6.2: Penalty (at BER of 107?) for a random channel realization. 
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Figure 6.3: Penalty (at BER of 10^?) as a function of the ground reflection amplitude. 


6.2 Suboptimal Antenna Placements and XPD 


Installing the Tx- or Rx antennas with suboptimal separation will negatively affect the per- 
formance of the MIMO system since it lowers the condition number of the channel matrix. 
There is thus a penalty associated with 7 Z 1 (see (3.16)) and that is plotted for a flat channel 
(b = 0) in Figure 6.5(a) for a 2x2 SS-SP LOS MIMO as a function of 7. The system is not 
so sensitive to changes in ņ if it is close to one. For y = 0.9 a penalty of less than one dB 
is observed but the penalty grows quicker as 7 becomes smaller. This results is not in good 
agreement with Figure 3.14 since direct matrix inversion and SVD are fundamentally different 
approaches (see Appendix B for further clarification). 

In Figure 6.5(b) the penalty is given as a function of XPD for 3 different 7 for a 4x4 SS-DP 
LoS MIMO. As the XPD and are lowered the penalty increases as expected. 
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Penalty vs. Filter Length. Two parallel SISO Ch. 
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Figure 6.4: Penalty (at BER of 10^?) as a function of filter length for the SISO representation 
of channel A and B. Here the un-penalized SNR is considered to be 10.5 dB while it is 7.5 
dB for the MIMO systems (array gain). 
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Figure 6.5: Penalties (at BER of 107?) for suboptimal antenna placement 7 and low XPD 
values. 


Added penalty is also expected for the frequency selective channels A and B when lowering the 
XPD. This added penalty caused by low XPD values should as well grow similarly independent 
of the channel realization. This can be seen in Figure 6.6 where the penalty is plotted against 
XPD for channels A and B. The penalty difference between the channels is seen to be almost 
constant at 5 dB. 
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Penalty vs. XPD. 4x4 SS-DP LoS MIMO 
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Figure 6.6: Penalty as a function of XPD for channel A and B. 


6.3 Decision Directed Algorithms 


It is important that an adaptive algorithm responds sufficiently fast to changes, that is that 
it has fast enough convergence speed. For WRLS and LMS it is apparent that the weighting 
factors and step-sizes will play an important role in determining the adaptation speed. Se- 
lecting a static channel and setting the SNR to the very high value of 95 dB the convergence 
speed can be seen by counting how many symbols it takes for the algorithm to reach its 
residual MSE floor. LMS and WRLS are analyzed in DD-mode, both started with the same 
channel estimate obtained using LS and training. 

Figure 6.7(a) shows the DD-LMS algorithm converging to its residual MSE of channel A. As 
the step-size is increased the convergence becomes faster, in fact if the step-size is doubled only 
half the symbols are needed to reach the noise floor. A trade off is that the actual MSE floor 
reached is higher if the step-size is larger. If the step-size is chosen too large the algorithm 
becomes unstable and might diverge. The LMS with the step-size of ums = 6 - 1074 needs 
about 25 thousand symbols to reach the MSE-floor. 

Figure 6.7(b) shows how the WRLS algorithm behaves for the same setup. Three different 
weighting factors are chosen, Aj = 0.98, A9 = 0.99 and Az = 0.995. A lower weighting 
factor means that the algorithm is able to adapt faster but with the penalty of a less refined 
autocorrelation matrix which means that the residual MSE reached will lie higher. A higher 
weighting factor will give a slower convergence but a more precise autocorrelation matrix 
resulting in a better convergence. This can be seen by looking at Figure 6.7(b). A3 results 
in slower convergence than using A; but a lower MSE floor. About 500 symbols were needed 
to reach the MSE floor for WRLS with A1. A huge improvement compared to the LMS 
algorithm. 


6.4 Blind Algorithms 


The convergence of blind algorithms such as CMA, MCMA or SCMA is quite sensitive to 
their step-sizes. The step-size has to be carefully chosen in order to avoid divergence of the 


85 


LMS and Ch, for Rx, 


T 
Si 
m 
[7 
0.5 1 1.5 2 2.5 3 
Symbol [105] 
LMS and Ch, for Rx, 
-40 
-50 
S -60 
Ww 
® 70 
-80 
0.5 1 1.5 2 2.5 3 
Symbol [10°] 
(a) LMS 
WRLS and Ch, for Rx, 
-50- 
—60 L 
Oo 
2 
u -70r 
-80r 
0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
Symbol [10°] 
WRLS and Ch, for Rx, 
m 
2, 
m 
[7 


05 1 15 2 25 3 35 4 #45 5 
Symbol [10°] 


(b) WRLS 


Figure 6.7: Convergence of the WRLS and LMS to their MSE-floor. 


algorithms. There is a trade off between convergence speed and the residual MSE floor, a 
large step-size value leads to a fast convergence speed but relatively high residual MSE or 
even instability and vice versa. In order to make a fair comparison, the same step-size is set 
for all the blind algorithms when comparing the performance. This step-size is chosen by 
simply searching for a suitable step-size for all the algorithms. A step-size of 5 - 109 is set 


for the following simulations. 
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The convergence of blind algorithms is simulated using a static channel (no mast swing) for a 
few ground reflection amplitudes (b factor). A larger b factor leads to a stronger inter-symbol- 
interference. All the blind algorithms converge for both Rx; and Rx? when b = 0 as shown 
in Figure 6.8(a), both channels behave exactly the same and the convergence performances 
are the same. However, when b = 0.5 convergence is not always reached. Comparing the 
performance for channel A and B (Figure 6.8(b) and 6.8(c) respectively) it can be seen that 
it takes a longer time for algorithms in channel A to converge but a lower residual MSE floor 
is achieved. Channel B converges faster for MCMA and SCMA but CMA does not converge 
for Rx». 

A much stronger Inter-symbol-Interference is introduced in Figures 6.8(d) and 6.8(e) by set- 
ting b = 0.9. Both channel A and B become too difficult for the STE to converge, convergence 
fails for such a strong ground reflection value and inter-symbol-interference. 


6.5 Mast Swing 


Mast swing will cause the channel to vary as the antennas move. A changing channel will 
require adaptive algorithms that can adjust fast enough to compensate in real time. At high 
carrier frequencies only a small change in the distance between the Tx and Rx will rotate 
the constellation. This rotation in the LoS paths will cause a fast change in the frequency 
selective channels. The mast swing velocities are given in the unit A/s so that the results are 
independent of the carrier frequency. Another reason for using that unit is that the affects 
of the velocity of the mast swing depends on how the phases between the ground reflection 
and direct paths are changing in relation to the actual mast movement. A change of one 
wavelength can be considered to be a change of 27 for those phases. Assuming a varying 
direct path but a static ground reflection then for the carrier frequency selected here (17 
GHz) the relation 1 A/s zz 0.02 m/s is valid!. 


6.5.1 Decision Directed 


Decision directed algorithms will diverge if they are fed with sufficiently many wrong symbol 
decisions. If such an algorithm loses track once, the errors will saturate and the algorithm 
will diverge to some false realization. If the algorithm adapts too fast then fewer decision 
errors are needed to get it to diverge, if the algorithm on the other hand adapts too slowly 
more decision errors will be made in a varying channel resulting again in divergence. Thus 
for a particular system setup and channel there is an optimum selection of parameters for the 
algorithm chosen. As the velocity of the mast increases it gets harder to adapt sufficiently 
fast until the algorithm diverges for any parameter selection. To evaluate the performance of 
the algorithms the squared-errors (SE) are plotted against the number of symbols received. 
A moving average filter is applied to the SE to smooth out the symbol-to-symbol variations 
making it a moving mean-squared-error. The initial channel estimate is obtained by sending 
training symbols in the beginning. For a BER of 107? a MSE of about -7 dB is needed for 
16QAM and a pure Gaussian noise channel. 


The antennas move with slightly different velocities because the mast is tilting. 
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WRLS 


The DD-WRLS algorithm is applied for a 2x2 SS-SP LoS MIMO system for different mast 
swing velocities. Three weighting factors were used to get a comparison of the performance. 
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Figures 6.9(a) and 6.9(b) plot the evolution of the SE for channel A and B respectively for 
both receivers. The velocity of the mast is about Vms = 182 A/s which means that the 
constellation is rotated seven times during the entire simulation. That produces about seven 
peaks which vary a little bit in shape due to the varying tilt of the masts. These peaks are 
introduced when for example a notch moves into or out of the bandwidth of the channel 
forcing an abrupt change in the equalizers. 

There are three weights A = 0.98, A = 0.99 and A = 0.995 chosen, and they are seen to give 
different results. The lower the weight the faster the algorithm will adapt (smaller exponential 
window), thus it can become too susceptible to erroneous decisions. This is seen to happen 
for the lower weights which loose track when there are decision errors (high MSE). 
Choosing the weighting factor too large will make the algorithm too slow and at high velocities 
it will lose track as can be seen in Figures 6.10(a) and 6.10(b) where A = 0.999 performs 
notably worse. None of the weights are able to keep track of the channel when the mast 
velocity reaches Vms = 286 A/s for channel B. 


LMS 


The LMS algorithm is much more sensitive to the velocity of the mast since its convergence 
speed is much slower. Figures 6.11(a) and 6.11(b) show the SE for channel A at velocities 
Vus = 4 A/s and Vms = 8 A/s respectively. A training sequence is transmitted every 
10° symbols. A new channel estimate for every training sequence is obtained using the LS 
algorithm. The step-size of ms = 6- 1074 is seen to be unstable since it diverges completely 
even after receiving a training sequence”. The step-size of ums = 5- 107? is too slow even 
for this moderate mast swing velocity of Vms = 4 A/s. At a velocity of Vms = 8 A/s all the 
step-sizes diverges for Rxg. As can be seen the adaptation speed is greatly inferior to the 
WRLS algorithm. The step-size of iius = 0.0006 is seen to be unstable while the lowest 
step-sizes are too slow even for a moderate mast swing. 


Higher Order Modulation 


A larger constellation requires a more precise channel estimate and thus more stringent re- 
quirements on the DD-algorithm that is compensating for mast swing. Thus a reduced tol- 
erance to mast swing is expected for 256QAM compared to 16QAM. Figures 6.12(a) and 
6.12(b) show the evolution of the SE for the WRLS algorithm. It fails to track somewhere 
between Vms = 46 and Vms = 91 A/s which are much lower velocities than for 16QAM. 


SISO Systems 


The SISO channel is much easier to track since the channel conditions can not become as 
complex as in the MIMO case. Figure 6.13(a) shows how the LMS adapts to mast swing for 
five different step-sizes. The two largest of them are seen to be unstable and the two smallest 
do not have sufficient adaptation speed. Only one step-size of ms = 107? is able to track 
at the mast swing velocity of Vms = 46 A/s. The evolution of the SE at Vms = 2.3 kA/s 
is plotted for three different weights of the WRLS algorithm in Figure 6.13(b). The weight 
A = 0.98 is able to track at this velocity. 


?Only every 1:1000 symbol is plotted. 
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Figure 6.9: Evolution of the squared error (SE) for three different weights in the WRLS 
algorithm in DD-mode in the presence of Mast Swing. Vms = 182 [A/s]. 


It is clear that mast swing is a much easier to compensate for in SISO systems compared to 
MIMO. 
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Figure 6.10: Channel B: Evolution of the squared error (SE) for three different weights for 
the WRLS algorithm in DD-mode in the presence of Mast Swing. 


91 


LMS Rx, and Ch, - Vus" [Vs] 


Symbol Nr. [10°] 
LMS Rx, and Ch, - Vus" [A/s] 


0 5 10 15 20 
Symbol Nr. [10°] 


(a) Channel A 


LMS Rx, and Ch, - Vyg=7.6 [s] 


Symbol Nr. [10°] 
LMS Rx, and Ch, - V,,,=7.6 [s] 


SE [dB] 


Symbol Nr. [10°] 


(b) Channel B 


Figure 6.11: Channel A: LMS with a training sequence after every 107% Tx symbols. 
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Figure 6.12: Channel A: Evolution of the squared error (SE) for three different weights for 
the WRLS algorithm in DD-mode in the presence of Mast Swing. 
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Figure 6.13: Evolution of the squared error (SE) for a SISO representation of channel A. 
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6.5.2 Blind Algorithms 
Flat channel 


For a flat channel, mast swing behaves as a phase rotator that constantly rotates the phase 
of received signals, this results in a constant carrier frequency offset between the transmitters 
and receivers. Figure 6.14(a) illustrates the signal constellation of the equalized signals. The 
plots to the left illustrate the constellation of signals along with the soft decisions of the 
STE without suffering any mast swing. The plots to the right illustrate the impact of carrier 
frequency offset that is caused by mast swing. 


The rotation shown in Figure 6.14(a) is caused by a carrier frequency offset and it is usually 
compensated for by using a carrier frequency recovery loop|31]. However, in Figure 6.14(b) 
and Figure 6.14(c) the algorithms MCMA and SCMA are able to accomplish blind equaliza- 
tion and carrier phase recovery simultaneously without the use of a carrier frequency recovery 
loop?. 


The convergence of three blind algorithms is compared in Figure 6.5.2. In order to compare 
the convergence of CMA with other blind algorithms, phase recovery is used for CMA, all 
three algorithms use the same step-size. It can be seen that all algorithms have converged 
after 20,000 symbols and that the residual MSE floor is stable even when the mast swing 
velocity equals to Vms = 131 A/s. The conclusion is that the convergence of blind algorithms 
is not affected by mast swing if there is a flat channel. 


Frequency Selective Channel 


Mast swing varies the MIMO channel so that the notch placements moves inside and outside 
the bandwidth of the signals. After convergence, blind algorithms can track the channel for a 
slow varying frequency selective channel but with a relatively high residual MSE floor. One 
solution is that the system uses blind algorithms to initialize the channel estimation and 
then it switches into DD (Decision Directed) mode to track the faster channel changes and 
improve the residual MSE simultaneously. Figure 6.16 shows the convergence of different 
blind algorithms with a mast swing velocity of Vms = 131 A/s and the ground reflection 
factor b = 0.9. 

As seen from Figure 6.16, all algorithms converge and the residual MSE floor is relatively 
stable for certain ranges (e.g. 100,000 to 300,000 symbols). The blind algorithms are able to 
track the channel to some extent but often with a relatively high MSE. The blind algorithms 
never diverge or lose track since they do not rely on previous decisions. They can thus over- 
ride bad channel conditions and then converge again when the channel conditions become 
better. Compared to LMS or WRLS, blind algorithm have a very high MSE which can result 
in decision errors. Thus it is desirable to switch to DD (Decision Directed) mode to improve 
the channel tracking and MSE as shown in Figure 6.5.1 for example. 


Figures 6.17(a) and 6.17(b) shows square error curves when the step-size is set to jj = 1-1075 
for channels A and B respectively. CMA and SCMA can become unstable and fail. MCMA 


If no other carrier frequency offset is introduced. 


95 


CMA: Vus = 131 Us 


Quadrature 
Quadrature 


0 0 
In-Phase In-Phase 


(a) Equalized signals using CMA 


MCMA: Vus = 0 Ms MCMA: Vus = 131 Ms 
4 4 
. a ds t4 é E EM 
+ $4 s o* € et 
2r 2r 


Quadrature 
o 

Quadrature 
o 


-2r -2r 
kirE] < $143 
—4 1 -4 ñ 
-5 0 5 -5 0 5 

In-Phase In-Phase 

(b) Equalized signals using MCMA 

SCMA: Vus -0JXs SCMA: Vus = 131 Us 

4 4 
Seeks | | REE 
2 | 2 ES 


D&$3) 1. FREY 


Quadrature 
o 

Quadrature 
o 


++ Y o d Os PHY 


-2} 1 - 

sët 3 SHH 

4 0 5 pi 0 5 
In-Phase In-Phase 


(c) Equalized signals using SCMA 


Figure 6.14: The compensation of carrier frequency offset using different blind algorithms 
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Figure 6.15: Convergence for a time varying flat channel 


is though still able to converge for this step-size. This goes to show that the algorithms are 


very sensitive to the step-size. 
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Figure 6.16: Square error curve for a time varying frequency selective channel. 
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Figure 6.17: Square error curves for a time varying frequency selective channel with a larger 
step-size. 
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6.6 Phase Noise 


The amount of phase noise power that a communication system can tolerate is highly depen- 
dent on the actual implementation of the PLL (or phase tracker). The power spectrum of the 
phase noise is also important to consider as well as its correlation between Rx branches in a 
MIMO system. This complicates simulations that aim to give a general result. Comparing 
the performance for different channel conditions using a fixed implementation can though be 
enlightening. 


The phase noise is simulated as 1/f? noise and it has a power of -66 dB[rad?]/Hz power at 
100 kHz from carrier (see Figure 5.11). 


6.6.1 Decision Directed 


Figures 6.18(a) and 6.18(b) show the MSE for different powers of the phase noise for 2x2 
SS-SP LoS MIMO and a corresponding SISO system respectively. The parameter Wp is 
added to the spectrum shown in Figure 5.114. The simulations for both the MIMO and SISO 
use the same implementation of STE and phase tracking and the phase noise is perfectly 
correlated between the receivers. 

The phase tracker should compensate for the phase noise before the STE so little noise 
enhancement is observed for different b factors other than the usual white Gaussian noise 
enhancement. The step-size is seen to play only a minor role due to the fact that the channel 
is static. 

As can be seen from Figure 6.6.2 the MIMO and SISO systems behave similar for fully 
correlated phase noise. For a differential phase noise the SISO system has the same tolerances 
but a MIMO system will perform worse for the setup illustrated in Figure 5.1. Applying 
differential phase tracking, that is a separate phase rotator in front of every equalizer, can 
compensate for this. Differential phase noise is applied to the received signal in Figures 
6.19(a) and 6.19(b) with and without a phase tracking respectively. Since the phase tracking 
is not designed for differential phase noise the system is seen to be rather sensitive. For 
very low phase noise power the simple phase tracking performs marginally better than no 
phase tracking but is seen to saturate the error for higher noise powers. Now the channel will 
enhance the phase noise power since the phase tracker does not compensate in front of the 
inter-stream-cancelers and thus the performance is sensitive to the b factor. 


6.6.2 Blind Algorithms 


In Figures 6.20(a) and 6.20(b) the convergence and residual MSE for the blind algorithms 
can be seen at different phase noise power levels for fully correlated phase noise. Behind the 
STE is a phase rotator that rotates the symbols back using the information from the last 10 
symbols, this is done to obtain a SE curve. At high phase noise powers the phase will actually 
vary quite a bit over 10 symbols raising the achievable residual MSE. 

CMA is unaffected by the phase noise while MCMA and SCMA are somewhat affected by it® 


“Won is added equally to all frequencies. 
5MCMA and SCMA both use not only amplitude but some phase information as well and are thus not 
completely unaffected by phase noise. 
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Figure 6.18: Channel A: The MSE for different power of the phase noise spectrum for fully 
correlated phase noise. 


which can be seen by looking at Figures 6.20(a) and 6.20(b). Here MCMA and SCMA have 
a higher residual MSE than CMA does. 

Figures 6.21(a) and 6.21(b) show the response of the blind algorithms to differential phase 
noise. Since the phase of the inter-stream-interference is constantly changing the STE can not 
converge fast enough at high phase noise powers. The blind algorithms are able to converge 
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Figure 6.19: Channel A: The MSE for different power of the phase noise spectrum for differ- 
ential phase noise. 


for Wo, = —5 dB but fail for W,, = 0 dB for this specified channel and step-size. 
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Figure 6.20: Convergence of the blind algorithms in the presence of fully correlated phase 
noise. 
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Figure 6.21: Channel A: Evolution of the SE in the presence of differential phase noise for 
the blind algorithms. 
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Chapter 7 


Conclusions and Discussions 


The following conclusions have been drawn from the numerical analysis: 


Higher penalties can be seen for frequency selective LoS MIMO channels compared to 
SISO. The channels also place a higher demand on the equalizers, both their struc- 
ture and their size. Equalizing a LoS MIMO channel is thus harder than its SISO 
counterpart. 


For a 2x2 SS-SP LoS MIMO system the penalty is below 1 dB if the separation between 
the antennas at the Rx and Tx are within 10% of the optimal. 


Comparing a 2x2 SS-SP LoS MIMO to a SISO system with a similar setup, tolerable 
mast swing velocity can be at least ten times lower for the MIMO system. 


For the equalizers WRLS outperforms LMS greatly for a time varying channel. The 
mast swing velocity that can be tolerated is about 100 fold greater for WRLS than it 
is for LMS. 


For the channels and system setup studied in this thesis the maximum mast swing 
velocity that can be tolerated lies around 5 m/s for WRLS. 


Blind algorithms may not convergence in the presence of strong ground reflection. How- 
ever, since mast swing introduces a time varying channel it can assist the blind algo- 
rithms by moving them out of an unfavorable channel condition. 


In the absence of phase tracking, convergence is difficult to achieve in the presence of 
differential phase noise for the blind algorithms. Having decision directed phase tracking 
the blind algorithms need to be able to converge before phase tracking can be initiated. 


MCMA and SCMA are able to compensate for carrier frequency offset in the kHz range. 


CMA is more tolerant to phase noise than MCMA and SCMA who use some phase 
information. 


SCMA has lower computational complexity without loss in performances. 


Blind algorithms may converge to the same source signals. This issue can be avoid by 
adding a cross-correlation term into the error functions. 
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Implementing SS LoS MIMO will put higher requirements on the structure of the system than 
on a SISO system, especially on signal processing. Channel conditions have a higher degree 
of freedom which means that it can be harder to tolerate unfavorable channel conditions. For 
wireless microwave links that need to be almost always operational this translates to larger 
system margins with some associated cost. 

Blind equalization is seen to be less reliable for MIMO than SISO and especially in the 
presence of differential phase noise and strong ground reflection. This might become a possible 
issue for future systems. 


7.1 Further Studies 


There is much room for further studies in this field. There are many components that are 
needed to build even a basic communication system. Their interaction as well as their im- 
plementation is challenging to handle. In order to be able to interpret and focus on specific 
phenomena in this thesis some simplifications were made to relax the requirement on the 
design. These simplifications should be dropped for future research. 

By dropping the approximation of a synchronous system the requirements on the system are 
increased. Carrier frequency offset recovery, symbol timing estimation, phase adjustment and 
tracking between different parts of the system and frame synchronization would for example 
be needed. These are all candidates for further studies and simulations. 

The simulations here have shown that in the presence of differential phase noise there can be 
convergence issues with blind equalization if there is no phase tracking available. This could 
represent an interesting research topic. 

Blind algorithm may converge in order to recover the same source signals at many outputs. 
Researches for improving and assuring the output separation are needed. 

Symbol timing estimates, fractionally spaced equalizers, non-linear equalizers and coded sys- 
tems are a few obvious extensions to the work already done. Another thing that might be 
improved on is the channel modeling. With some real measurement data for LoS MIMO 
channels some statistical representation of the channel could be made which should allow for 
more typical results. 


So to summarize the future work suggested is as follows. 
e Handling differential phase noise 
e Improving the convergence of the blind algorithms 
e Fractionally spaced equalization and non-linear equalization 
e Symbol timing estimates 


e Using real channel measurements when available 
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Appendix A 


Derivations 


A.1 Singular Values For 2x2 LoS MIMO 


The singular values for a MIMO system are obtained by calculating the eigenvalues of W 
where 


uas HH" if N<M 
«| HEH if N>M 
For a 2x2 LoS MIMO we have 


Hs = | exp(jkd11) exp(jkdi2) | E 
á exp(jkd21) exp(jkd22) 
u H 2 ejk(d11—d21) + edk(d12—d22) 
W = Hros Hios m | ejk(d21—d11) af ejk(d22—d12) 2 


The eigenvalues are found by solving the following equation 
det (W — cI) 2 0 
This yields 
(2 = c? E gi (dii +d22—d12—d21) = egi k(di2+da1—d22—do2) zx 


which gives 
c=2>+|[2+ 2cos(kA)] 2 (A.1) 


where A = d11 + do? — dig — d21. The singular values are now the square roots of the ordered 
eigenvalues c. The condition for optimality for the 2x2 MIMO system is that dj, — d12 = 3, 
thus we have Aopt = £z. For 7 defined in (3.16) the eigenvalues in (A.1) become 


e 


c= 2+ [2 + 2cos(rn? 


Thus the singular values can be expressed as 
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1/2 
gie E + [2+2 cos(77”)| ui 


1/2 
09 = E — [2 + 2cos(zi?)] i (A.2) 
And finally the condition number is 
A 
is 


A.2 Derivation of the normal equations 


For simpleness, we only proof the least square based on the output of the first receiving 
antenna. It can be extended into any number of receiving antennas. We use y; to express 
the output of the equalizer at the first receiving antenna and we are assuming there are N 
samples of the received signal. 


M-L 


ilk) =D riw; (A.3) 


j=1 
For all k = 1,2,...,N. M, L is the number of Rx antennas and equalizer filter length re- 
spectively. rk j, w; is the k, jth element in R matrix and jth element in column vector Wj. 


ew, (k) = &(k) — m (k = 1,2,...,.N). (A.4) 
j=l 


So to minimize the function under LS algorithm. 


N 
¿=D dy (h). (A5) 
kel 


where N is the number of samples from the first receiving antenna output before the equal- 
ization operation. e is minimized when its gradient with respect to each element of vector 
Wi is equal to zero. 


N 
Oe Oew, (k) 
=? AS y ] 
since "T 
ew, (k) = s1(k) - So rg;w; (k=1,2,...,N) 
j=l 
the derivative is 
dem M-L 
OW, i 2 Th 
j=l 


Substitution of the expression for the residuals and derivatives into the gradient equations 


obtains: 
N M-L M-L 
OE 
= -2 > ò rk j(sı(k) > Tg us) =0 
OW; 
=] j=l i=l 
After re-arrangement, the normal equation is re-written as: 
N MLN M-L N 
> , Tk j' Tki Wi = , ò Tk j’ sı(k). 
j=1 i=1 k=1 j=1 k=1 


The normal equation is written in matrix notation as: 
(RTRW,- R” S; 


Thus (4.16) is proved. 
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(A.T) 


(A.8) 


Appendix B 


Singular Value Decomposition 


For every N x M matrix H one can write 


H = UZV" (B.1) 


where U and V are M x M and N x N unitary matrices respectively. X is a M x N diagonal 
matrix of singular values o; of H. These singular values have the property that ok = VAk 
where A, is the ith largest eigenvalue of HH". If the transmission system is designed in 
such a way that it has transmitter precoding an receiver shaping that is if the sequence to be 
transmitted is multiplied with UP and the received sequence with V the following result is 
obtained [20, ch. 10] (see figure 3.2 for reference) 


r = (sH+n)V 

(sUXVH + n)V 

(x UPUXVH + n)V 

= xU'"UXVPv + nV 

= x% +ñ (B.2) 


where ñ = nV has the same distribution as n since V is a unitary matrix (VĦV = Im). The 
derivation above shows that with the proper precoding and receiver shaping a MIMO channel 
can be parallelized into Ry SISO channels which do not interfere with each other and are 
only connected through the total power constraint. 


B.1 Example: 2x2 MIMO 
Using SVD the SNRs for a 2x2 MIMO are given by (seen from (B.2)) 
o? oz 
SNR; = oi SNR; = oi (B.3) 


where o? is the noise power. Equalizing using direct matrix inversion will though give a 
different result as is shown here 
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r = sH+n 

y = rH ! -2 HH! ; nH^! 
= s+nVy U! 
= s4üX tU! 


Now the SNRs can be written as 


SNR 1 o? ae 
1 — pe: 
02101112/0% + o2|Uai|?/05 o2 03/U11l? + of |Uar|? 
2 2 
O O 
SNR» 2 l 


0? c2|Un + 02|Us1?? 


The SVD and direct matrix inversion approaches are thus not comparable. 
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